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FOREWORD 


This is a progress report on the research project, "Analysis and 
Computation of Internal Flow- Field in a Scramjet Engine," for the period 
ended December 31, 1989. Special attention during this period was directed 
to "Effects of Nose Bluntness and Shock- Shock Interactions on Blunt Bodies in 
Viscous Hypersonic Flows . " The work was supported by the NASA Langley 
Research Center (Computational Methods Branch of the Fluid Mechanics 
Division) through the grant NAG-1-423. The grant was monitored by Drs. Ajay 
Kumar and J. Philip Drummond, FLDMD- Computational Methods Branch, Mail Stop 
156. 
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EFFECTS OF NOSE BLUNTNESS AND SHOCK- SHOCK INTERACTIONS 
ON BLUNT BODIES IN VISCOUS HYPERSONIC FI0WS 

By 

D. J. Singh 1 and S. N. Tiwari 2 

SUMMARY 

A numerical study has been conducted to investigate the effects of blunt 
leading edges on the viscous flow field around a hypersonic vehicle such as 
the proposed National Aero-Space Plane. Attention is focused on two specific 
regions of the flow field. In the first region, effects of nose bluntness on 
the forebody flow field are investigated. The second region of the flow 
considered is around the leading edges of the scramjet inlet. In this 
region, the interaction of the forebody shock with the shock produced by the 
blunt leading edges of the inlet compression surfaces is analyzed. Analysis 
of these flow regions is required to accurately predict the overall flow 
field as well as to get necessary information on localized zones of high 
pressure and intense heating. The results for the forebody flow field are 
discussed first followed by the results for the shock interaction in the 
inlet leading edge region. 

The forebody is modeled by slender cones and ogives with a spherically 
blunted nose. A combination of Navier-Stokes and parabolized Navier-Stokes 
equations is used to compute the flow field. The influence of entropy layer 
thickness on the extent of the leading edge effects is also considered. The 
extent of downstream effects of leading edge thickness are determined at Mach 


Graduate Research Assistant, Department of Mechanical Engineering and 
Mechanics, Old Dominion University, Norfolk, Virginia 23529. 

^Eminent Professor, Department of Mechanical Engineering and Mechanics, Old 
Dominion University, Norfolk, Virginia 23529. 


iv 


numbers of 10 and 20 for cone angles of 5° , 10°, and 20°. Three values of 
nose bluntness are considered with the smallest nose blunting (0.0025m) 
representing the sharp cone/ogive. Depending upon the flow conditions and 
the geometry, significant differences have been found between the sharp and 
the blunted bodies; the differences persist as far as 236 nose radii 
downstream. Also, the bluntness effects decrease with increasing cone angle. 
The results show that the wall quantities are not affected much by the 
inclusion of high temperature effects through equilibrium chemistry. 

For the flow region around the inlet, depending upon the Mach number and 
the angle of attack, the forebody shock can interact either with the blunt 
cowl leading edge shock or with the shock produced by the blunt leading edges 
of the swept sidewall compression surfaces of the inlet. For the interaction 
at the cowl leading edge, the forebody shock is assumed planar and the cowl 
is modeled by a two-dimensional cylindrically blunted wedge of infinite 
width. Use of the full Navier- Stokes equations is made on the cowl forebody 
and the thin- layer Navier- Stokes equations are suitably modified for space 
marching on the cowl afterbody. The results of the study show that the flow 
around the cowl is significantly altered by the impinging shock. The peak 
value of pressure is found to be nine times and heating rates eight times the 
stagnation point value for the unimpinged case at Mach 8.03. The peak values 
were slightly lower for Mach 5.94 calculations. The use of solution adaptive 
grids is made to properly resolve the flow field. The peak heating rates for 
the unadapted grid was four times the unimpinged stagnation point value 
versus eight times with the adapted grid. A preliminary study was also 
conducted to determine the influence of shock- shock interaction on the blunt 
leading edges of the swept sidewalls of the inlet. For this configuration, 
the flow field is fully three-dimensional. A three-dimensional thin- layer 
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Navier - Stokes code had been used to calculate the flow field. The peak 
pressure for this case is found to be 2.25 times and the peak heating three 
times the unimpinged stagnation values. The results of the study are 
compared with the available experimental and numerical results. This study 
presents the first full three-dimensional analysis of the shock-shock 
interaction on the swept inlet sidewalls and the accurate numerical solution 
to date using full Navier-Stokes equations for the two-dimensional 
interaction on the cowl leading edge. 
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Chapter 1 

INTRODUCTION 


There is a renewed interest in the hypersonic flow regime after a gap of 
over a decade primarily because of the proposed National Aero-Space Plane. This 
hypersonic aircraft is to be powered by an airbreathing engine which should be 
very closely integrated with the airframe to avoid severe drag penalties. The whole 
forebody of the vehicle is used to precompress the airflow before it enters the engine 
inlet. The air is further compressed inside the inlets and the supersonic combustion 
takes place inside the scramjet combustor. The burned gases are expanded through 
the nozzles as well as on the aft portion of the undersurface. Integrating the air- 
frame with the propulsion system for the hypersonic aircraft makes the traditional 
wind tunnel guided design not only difficult but very expensive. However, due to 
recent advances in computer architecture and efficient algorithm, it is now possi- 
ble to solve more and more complex problems numerically. The Computational 
Fluid Dynamics (CFD) has emerged as an extremely valuable engineering tool in 
aerodynamic analysis and design. A wide variety of very complicated flows are be- 
ing simulated by computers. At the NASA Langley Research Center, a great deal 
of numerical and experimental efforts are directed toward a better understanding 
of the complex flow field in different regions of the Aero-Space Plane. Numerical 
modeling of the flow field has proven to be a valuable tool for getting better insight 
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into the complex nature of these flows. This approach in conjunction with ongoing 
experimental program provides an effective analysis capability. 

The forebody geometry is relatively simple and fixed and can be modeled 
by long slender bodies of revolution with blunted nose. A reusable vehicle will 
probably not have an ablative surface. Thus, the forebody flow field analysis appears 
fairly straight forward. But, in addition to the prediction of wall properties, such 
as skin friction, heat transfer and pressure, an accurate estimate must be made of 
the entire flow field in order to predict the mass and the momentum flux entering 
the inlet. The thickness of the boundary layer at the inlet face is important for 
predicting the inlet performance. It is, therefore, necessary to incorporate realistic 
geometrical modeling and flow conditions in the analysis of the forebody. One of 
the geometrical features which can have significant effect on the entire flow field is 
the nose bluntness. In hypersonic flows over slender bodies, the influence of nose 
bluntness can be significant for hundreds of nose dimensions downstream. Thus, it 
is important to develop efficient and reliable methods for predicting the effects of 
blunt nose on the entire flow field. 

In most studies involving flow past slender bodies, it is assumed that the 
leading edge is infinitely sharp. In practice, it is impossible to manufacture a body 
with zero leading edge thickness as well as all bodies must have some finite bluntness 
in order to sustain the heat generated at hypersonic speeds. The leading edge heat 
transfer can be decreased by increasing the leading edge thickness. However, this 
results in a drastic increase in drag. With the advancement in material technology, 
it is possible to reduce the leading edge thickness without encountering intolerable 
heating rates in the vicinity of the nose. But the extent to which this bluntness 
influences the development of flow downstream is not known precisely for hyper- 
sonic flows. In general, such a calculation is very difficult because of the various 
competing effects such as pressure interaction, vorticity interaction, and curvature 
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parameters. To calculate the flow over blunted slender bodies, two sets of equations 
are utilized in the present investigation. For the nose section, the use of full Navier- 
Stokes equations is necessary due to the presence of the subsonic region; but slightly 
downstream of the nose region on the body, the flow becomes predominantly super- 
sonic and the space marching parabolized Navier-Stokes equations can be used. A 
shock fitting algorithm for both the Navier-Stokes and parabolized Navier-Stokes 
equation is used. The real gas effects which might be present at hypersonic speeds 
are also considered. 

Another impact that a blunt leading edge can have is on the inlet flow 
field. The optimum propulsion efficiency occurs when the forebody shock impinges 
on the engine cowl lip so that all the precompressed air is captured by the inlet. 
Depending upon the Mach number and the angle of attack the forebody shock can 
interact with the shock from the compression surface leading edges of the inlet. 
This interaction results in a very complex flow field and can result in large increase 
in pressure and heat transfer over certain localized regions. The large temperature 
gradients cause thermal stresses which could result in structural failure. In order to 
relieve the influence of thermal stresses, some form of active cooling is needed. To 
determine the cooling requirements, pressure and heating rates on the body need 
to be predicted accurately. Such phenomena also occur on space shuttle and other 
high speed vehicles at the point where the bow shock from the nose interacts with 
the wing leading-edge shock. 

For accurate predictions of shock interference patterns and aerothermal 
loads, it is necessary to employ modem algorithms. The classical central difference 
numerical algorithms have less than desired resolution in the shock regions. In 
recent years, the upwind biased algorithms have become a popular alternative to 
the central difference schemes. These methods model the characteristic nature of the 
equation. The information at each grid point is obtained from the direction dictated 
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by the characteristic theory. The particular upwind method used here is the flux 
vector splitting due to van Leer which is used to solve the compressible Navier- 
Stokes equations with upwind biased differencing for the inviscid terms and central 
difference approximation for the viscous terms. The cowl forebody is modeled by a 
cylindrically blunted wedge of infinite width which allows a two-dimensional analysis 
of the flow field. An adaptive grid technique is used to properly resolve the flow 
field features and to accurately predict the surface heat transfer. The use of full 
Navier-Stokes equations is made on the cowl forebody and the thin layer Navier- 
Stokes equations are suitably modified for the space marching on the cowl afterbody. 
For the forebody shock interaction with the inlet sidewall, the compression surface 
leading edge is modeled by a swept wedge. In this case, the line of intersection of the 
forebody shock and the compression surface leading edge shock is curved making 

the flow field fully three-dimensional. 

The objective of the present investigation is to develop an accurate solu- 
tion procedure for determining the zone of influence of the leading edge thickness 
and to calculate the shock interference phenomenon with an accurate prediction of 
aerothermal loads on the blunt leading edges of slender bodies in viscous hypersonic 
flows. The high temperature effects which might be present at hypersonic speeds 
are accounted for by using an equilibrium gas model. The influence of entropy 
layer on the leading edge effects is also addressed. To accurately predict the shock 
interactions, the importance of a suitable grid is discussed. Two types of shock 
interference patterns are considered. The analyses ranges from a few simplified ap- 
proaches to the solution of the combination of full Navier-Stokes and parabolized 
Navier-Stokes equations. The solution procedures are validated by comparing the 
present results with the available experimental and numerical data. 


Chapter 2 

LITERATURE SURVEY 


In this chapter, the theoretical and analytical works carried out in 
the area of the effects of nose bluntness and shock-shock interactions are reviewed. 
The effect of blunt leading edge on flat plates and slender bodies is considered first 
and then the available work on shock-shock interactions is reviewed. 

Early studies on the effects of nose blunting are summarized in classical 
books on hypersonic flows [l] - [5]*. They deal with approximate approaches such as 
the Newtonian theory, small disturbance theory, and constant density assumptions. 
Some empirical correlations have been developed on the basis of experimental data. 
But these approaches have limited validity and are very restrictive in nature. 

Some specific studies of nose bluntness effects have been carried out for 
flat plates [6] - [10]. Hammitt and Bogdonoff (6] carried out experiments on blunted 
flat plates. A wide range of Mach number and Reynolds number were considered 
for various leading edge thicknesses. It was found that the flow over a flat plate 
with thick leading edge is essentially inviscid but viscosity plays an important role 
for thin leading edges. Also the effect of leading edge is felt for several thousand 
leading edge thicknesses downstream. Vas et al. [7] did an experimental study of 
the shock shape and surface pressure distribution about thick two-dimensional and 
axisymmetric bodies. The shock shape was found to be parabolic about flat plates 
‘The numbers in brackets indicate references. 
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with flat and cylindrical leading edge for < 10. The pressures on the cylindrical 
sections can be predicted reasonably well by the Newtonian theory and the pressure 
coefficient ratio on the surface was found to be higher for two-dimensional bodies 
than for axially symmetric bodies. Allegere and Bisch [8] studied the effects of the 
angle of attack and leading edge on the flow about flat plates at Mach 18. The 
Reynolds number considered were 70-15000. It was found that the wall pressure 
increases with increasing leading edge thickness and angle of attack. 

Cheng et al. [9] did a theoretical and experimental study of the leading 
edge bluntness and boundary layer displacement effects for the situation where both 
effects axe equally important. The interaction of the two effects was treated theoret- 
ically by extending the basic shock layer concept. In the outer inviscid flow, a model 
consisting of a detached shock wave and an entropy layer was introduced to account 
for the bluntness. In the boundary layer the approximate solution was found to be 
governed by the local flat plate similarity theory. Bradfield et al. [10] conducted 
an experimental investigation to ascertain the effects of leading edge bluntness on 
the laminar boundary layer of a flat plate. It was demonstrated that the effect of 
the leading edge bluntness is an important consideration in the formulation and 
development of the boundary layer on a flat plate at supersonic speeds. The ratio 
of leading edge thickness to the value of boundary layer thickness at the measuring 
stations may be taken as criterion of the magnitude of this effect. 

The flow field over slender bodies has been studied by various investiga- 
tors. However, most of these studies are limited to rarefied flows [11] - [14] and/or 
to a very small nose bluntness [15] - [28]. Vas and Sierchio [ll] performed an exper- 
imental study to determine the downstream effects of nose bluntness in hypersonic 
rarefied flows. The study was carried out at Mach 25 and freestream Reynolds 
number of 11000 per inch on 5° half angle sharp and blunted cones. Significant 
differences between the sharp and blunted bodies were noted on the surface and 
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flow quantities. McCroskey et al. [12] studied the leading edge flows over slender 
bodies in rarefied hypersonic flows. The primary attention was directed upstream 
of the merged layer. Based on the experimental data it was established that there 
is a competition between the displacement and merging effects in determining the 
surface pressure and heat flux. Feik et al. [13] conducted experimental studies on 
cones of 5°, 10° and 15° half angles at Mach 25 and Reynolds number 7000-15000 per 
inch. The effect of bluntness within the merged region and the flow downstream of 
the merged region was studied. Significant effects of the nose bluntness were found 
on the location of merging, the location of shock wave after the merging, and both 
the surface pressure and heat transfer within and considerably downstream of the 
merged region. 

Vogenitz et al. [14] conducted numerical experiments on rarefied flows of 
a monoatomic gas about slender cones and flat plates at Mach 10 and 25; cone 
half angle ranged from 3°-15°. The gas model was composed of hard sphere model 
which accommodates completely to solid surface and reflect diffusively. Leading 
edge flows on sharp slender cone was found to differ substantially from that on 
the flat plates. Burke and Curtis [15] measured static pressure distribution on the 
surface of a blunted 7.5 degree half angle cone for Mach 8-18 and Reynolds number 
of 1 x 10 6 per foot. It was found that the position of minimum pressure and the 
rate of pressure recovery depends upon cone angle. Bertram [16] examined the 
aerodynamic effect of nose blunting on a 10° cone, the blunting consisted of a plane 
cut normal to cone axis. The shock shape was found to have an inflection at about 
40-60 nose radii from the nose and the surface pressure reaches a minimum. For 
lower values of Reynolds number, the surface pressure data were not affected to 
any considerable extent by viscous effects. Dewey [17] studied the slender vehicles 
operating in high-altitude, high-temperature environment. The study showed that 
a regime of significant practical interest exists in which the effects of leading edge 
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bluntness and viscous interaction are of comparable magnitude; the lift and drag 
depend to a large extent on the relative effect of the leading edge bluntness and 
boundary layer displacement. 

Cheng and Pallone [18] investigated the inviscid leading edge effects in hy- 
personic flow on the basis of small perturbation theory. The results were obtained 
under the assumption of one dimensional adiabatic flow of an ideal gas and infinite 
shock strength. Various correlations for surface pressure and shock shape were de- 
veloped for several values of 7. Lees and Kubota [19] analyzed inviscid hypersonic 
flow over blunt nosed slender bodies. Based on similarity theory, correlations for 
shock shape and surface pressure for slender bodies were developed. It was found 
that the shape of the bow shock wave, flow field and wall quantities are dominated 
by the blunt nose over a downstream distance many times greater than the char- 
acteristic nose dimension. Cheng [20], on the basis of small perturbation theory, 
developed the laws of similitude for hypersonic, inviscid, real gas flow field over 
blunted flat plates, cones and wedges. The results for real gas flow over a cone at 
zero yaw agreed with the laws of similitude developed. Also the thickness of the 
singular region near the surface of a blunt region was estimated and its effects on the 
pressure prediction was determined. Lees [21] calculated laminar heat transfer over 
blunt nose bodies at hypersonic flight speeds for limiting cases of thermodynamic 
equilibrium and diffusion as controlling factors for the rate of heat transfer. De 
Jamette and Davis [22] developed a simplified method for calculating the laminar 
heat transfer over blunt bodies with only freestream Mach number and the ratio of 
specific heat as input parameters. The viscous problem was simplified by using the 
axisymmetric analog for three-dimensional boundary layers in conjunction with the 
Lees’ laminar heating rate relation. The results for 15° and 20° cones at freestream 
Mach number of 10.6, for various angles of attack, showed good agreement with 
experimental data. 
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In recent years, certain numerical and experimental efforts have been di- 
rected to investigate the effects of nose bluntness on slender bodies [23] - [28]. 
Stetson [23] conducted an experimental study on blunted cones to determine the 
leading edge effects on the boundary layer transition. Based upon the results it was 
concluded that transition depends upon the ‘entropy swallowing distance’ which 
in turn is a strong function of the leading edge thickness. Nowak [24] did an ex- 
perimental investigation at Mach 6.7 to determine the effects of Reynolds number, 
angle of attack and the nose shape on aerothermal environment of a 12.4° half angle 
cone. Three nose tip configurations were tested on the cone; a three inch radius, 
a one inch radius and a sharp tip. Heating rate distribution, surface pressure dis- 
tribution, shock shape and shock layer profiles were measured and compared with 
predictions. Surface pressure and heating rates when normalized by their stagnar 
tion values were found to be independent of the freestream Reynolds number. With 
the increasing bluntness, the surface pressure took longer distance to recover from 
the nose overexpansion. 

Zoby and Simmonds [25] did engineering calculations over hyperboloids, 
paraboloids, ellipsoids and sphere-cones. The heating and constant or variable 
entropy conditions were coupled with an approximate inviscid code. Results were 
obtained for various angles of attack and good comparison with experimental data 
was obtained. Thompson et al. [26] did a numerical study to assess the applicability 
of viscous shock layer code and various engineering codes for aerothermal predictions 
over slender spherically blunted cones. The numerical results were compared with 
experimental data from flight and ground based tests. The effects of nose blunting 
and angle of attack on drag and heat transfer were illustrated. Gupta et al. [27] 
obtained numerical solution for hypersonic laminar and turbulent flows over slender 
bodies. The results were obtained from the viscous shock layer equations at Mach 
8-15 for sphere-cones and hyperboloids. Detailed comparisons were made with 
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other numerical results and experimental data to assess the accuracy of viscous 
shock layer solutions. Zoby et al. [28] conducted a parametric study of laminar 
nonequilibrium heat transfer on slender vehicles. The study included the variation 
of altitude-velocity conditions, cone half angle and nose bluntness. The effect of 
wall catalycity to reduce the heating rates was also investigated. 

A great deal of work is currently underway on shock-shock interactions at 
various places due to the practical significance of the phenomenon in complex high 
speed flows. In relation to the National Aero-Space Plane, special attention is di- 
rected to study the shock-shock interaction at the cowl plate and the inlet sidewalls 
of the scramjet engine inlet. In order to avoid the structural failure, it is very crucial 
to calculate the thermal stresses arising in the scramjet inlet. To have an optimum 
design in terms of structural strength and cooling requirements, accurate calcula- 
tions of surface pressure and heating rates are of utmost importance. Due to the 
complexity of the problem, most of the earlier works are either experimental or of 
highly simplified model where empirical inputs have been used. The effects of Mach 
number, Reynolds number, impinging shock strength and sweep angle have been 
studied on surface pressure and heat transfer over flat plates, cylinders, hemispheres 
and various leading edge geometries. The experimental data has unusually large 
scatter due to difficulties in making accurate measurements over the narrow zone 
where the impingment occurs. Until recently, the numerical calculations were virtu- 
ally nonexistent. With the rapid advancement in computer hardware and improved 
algorithms, it is now possible to compute such complex flow field accurately. 

Teterin [29] conducted experimental study on the effect of incident shock 
on the flow field and force effects on the cylinder. The cylinder was mounted on a 25° 
wedge. In order to isolate the effects of shock-shock interaction from separation, 
the separation zone was eliminated by means of a fillet. Schlieren pictures were 
obtained and the pressure distribution was measured. Sukovstin and Shestova [30] 
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studied effects of bow shock impingment on a body located in the influence zone 
of primary body (creating the bow shock) . An approximate integral method was 
presented and results were obtained for Mach 3-20. 

Edney [31] in 1968 conducted a detailed study of the shock-shock inter- 
action phenomenon at Aeronautical Research Institute of Sweden. Three basic 
models, a 30 mm diameter hemisphere/cylinder, a 30 mm diameter flat faced cylin- 
der, and a 30° half angle cone/cylinder with 5 mm nose radius and 15 mm base 
radius were tested at Mach numbers of 4.6 and 7.0. Edney was the first to clas- 
sify these interactions into six basic patterns, which can occur when an extraneous 
shock interacts with the bow shock of a blunt body. The interactions could result in 
shock-boundary layer interaction (Type I, II, V), shear layers (Type III), supersonic 
jet (Type IV) or expansion fan-boundary layer interactions (Type VI). All of these 
interactions, except for Type VI, increase the surface pressure and heat transfer 
over a localized zone. The Type IV interaction is the most severe one, producing 
largest increase in pressure and heat transfer. The type VI interaction is the only 
one that causes the heat transfer and pressure to decrease due to production of an 
expansion fan. 

Keyes and Hains [32] continued the work of Edney for higher Mach num- 
bers, Reynolds numbers and various values of specific heat. All six types of interac- 
tion were studied. The planar shock waves were allowed to impinge upon the bow 
shock surrounding swept fins, hemisphere and wedges. The impinging shock wave 
was generated by a wedge. Heat transfer measurements were obtained using a phase 
change paint technique. Peak heating rates of 17 times and peak pressure of eight 
times the unimpinged stagnation point values were measured. The results showed 
that pressure and heating rates are strongly affected by Mach number, specific heat, 

impinging shock strength and geometry. 

Craig and Ortwerth [33] investigated the leading edge shock impingment 
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problem at Mach 15 in the Airforce Applied Physics Laboratory. This study was 
carried out to determine the pressure and temperature amplifications on the leading 
edge of hypersonic inlets. The peak heat transfer amplification on leading edge was 
found to be less than five. In this study the spacing between the pressure taps 
and the thin film thermometer was too large to give adequate resolution. The data 
obtained cannot be used due to lack of resolution and unavailability of information 
on the position and strength of the impinging shock. Wieting [34] conducted 
an experimental study on the shock wave interference over a cylindrical leading 
edge at Mach 6.3, 6.5 and 8.0. The model consisted of a three inch diameter 
cylinder. The results were presented for heat transfer and pressure distribution for 
a two-dimensional shock wave interference on a cylinder. This was considered to be 
representative of the cowl leading edge of the scramjet inlet. The local heat transfer 
rate and pressure were amplified up to ten times the unimpinged stagnation point 
values. 

Several semi-empirical approaches have been proposed to theoretically pre- 
dict the peak heating and pressure [35] - [40] but they rely on several empirical 
inputs which must be known apriori. Morris and Keyes [35] developed a computer 
program based on semi-empirical correlations to calculate the peak pressure and 
heating rates for all six types of interactions. The correlations require various in- 
puts such as width of impinging jet and length of transmitted shock. In general this 
information is not known apriori. Crawford [36] improved upon the Edney s graph- 
ical approach for predicating peak pressure and heat transfer by plotting a family 
of pressure deflection curves with the pressure ratio on a logarithmic scale and flow 
deflection on linear scale thus eliminating the need for iterations as Mach number is 
uniquely defined on this set of coordinates. Keyes and Morris [37] developed correla- 
tions for peak pressure and heating for shock/boundary layer interaction and shear 
layer attachments for two-and three-dimensional interactions. The Mach number 
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considered were 6 and 20 over a freestream Reynolds number range from 3.3-25.6 
million per meter and specific heat ratio of 1.4 and 1.67. The shock generator angle 
was varied from 5° - 25°. Bramlette [38] developed an approximate method for 
computing pressure for Type III and IV interference. The peak pressure obtained 
for various shock generator angles showed good comparison with the experimental 
data. Markarian [39] studied the problem of aerodynamic heating caused by shock 
wave and boundary layer interactions. It was found that heating rates for this 
type of interaction are significantly higher than normally predicted values; also the 
heating rates were found to be higher for laminar interaction than turbulent inter- 
action. Some empirical correlations were developed to predict the heating based 
on the inviscid pressure rise across the interaction. Bertin et al. [40] studied shock 
interference pattern for delta wing orbiters. The surface pressure and heat transfer 
distribution in the interaction perturbed region were obtained for velocities from 
1167 m/s to 7610 m/s. At high velocities, the Type IV interaction was found to 
exist for wing leading edge sweep as low as 27^. For this case, no locally severe 
heating rates were obtained. 

There are very few numerical studies available on the shock-shock inter- 
action due to the complexity of the problem. Tannehill et al. [41,42] presented the 
first successful numerical simulation of the interaction problem based on the explicit 
MacCormack method. White [43] solved this problem using the Reynolds-averaged 
Navier-Stokes equations (based on an implicit finite-volume method) and obtained 
results for Mach 6.0. Klopfer [44] conducted an extensive study for various Mach 
numbers using the second-order implicit TVD algorithm but the results for the 
peak heat transfer for the Type IV interaction were off by a factor of about 2.25, 
although the pressure was in close agreement with the experimental results. Perry 
et al. [45] used the Roe scheme to solve the Type IV interactions but did not give 
any heat transfer results. Moon et al. [46], in a recent paper, provided the results 
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for the Type III interference pattern. The study indicated that the flow for this 
type of interaction in the impingement region is turbulent rather than laminar. 

Most of the studies mentioned earlier to study the shock-shock interaction 
are two-dimensional in nature except [32], thus cannot be used for interaction on 
swept inlet sidewalls. The only related numerical study on three-dimensional inter- 
action was carried out by Holst et al. [47]. They solved the three-dimensional shock- 
shock interaction on an infinite cylinder using the explicit MacCormack scheme. 
The viscosity was assumed an order of magnitude higher than the actual viscosity 
to physically thicken the boundary layer and a very coarse grid was employed. The 
resulting solution was at best qualitative in nature. 

The literature survey clearly indicates the need for an additional extensive 
study in this important field of current interest. 


Chapter 3 


DOWNSTREAM EFFECTS OF 
NOSE BLUNTNESS 


In this chapter the effects of nose bluntness on the flow field and wall 
quantities over forebody of a generic hypersonic plane are considered. The calcu- 
lation of hypersonic viscous flow field past long slender axisymmetric blunt bodies 
is of prime interest to the designer of aerospace vehicles. A wide variety of flow 
conditions are encountered during the trans-atmospheric flight of these vehicles. As 
mentioned previously, the forebody of these vehicles must have finite thickness in 
order to sustain the pressure and heating rates encountered at hypersonic speeds. 
The forebody is modeled by spherically blunted cones and ogives. The physical 
model for hypersonic flow past a blunt-nosed slender body is shown in Fig. 3.1. 
The relative positions of boundary, entropy, and shock layers are indicated in the 
figure. The first physical problem considered is the flow past a blunt nosed slen- 
der cone. Numerical results have been obtained by employing a combination of a 
Navier-Stokes code and a parabolized Navier-Stokes code. These are discussed here 

briefly and some specific results are presented in Chap. 5. 

In the present study, 5°, 10°, and 20° half-angle cones are selected first with 
spherical nose tip as shown in Fig. 3.1. Three nose tips with radius of 0.0025m, 
0.025m and 0.05m are considered, with the smallest nose radius tip approximating 
the sharp tip. Blunting was accomplished by keeping the cone angle fixed and 
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increasing the nose radius. The bluntness effects are determined by comparing the 
solutions obtained from the blunted tips with that obtained from the sharp tip. As 
shown in Fig. 3.1, the origin of the coordinate system is at the virtual tip of the 
sharp cone with x axis along the symmetry line and y axis normal to it. Due to 
symmetry of the flow field at zero angle of attack, only half of the flow is calculated. 
The second geometry considered is a generic forebody (Fig. 3.2) with three nose 
tips of same radii as the cone. Here again the smallest nose radius approximates 

the sharp tip body. 

3.1 Entropy Layer 

Nose bluntness at hypersonic speed causes shock wave to be detached and 
curved in the nose region. The entropy of the flow increases across a shock wave, 
the stronger the shock wave, the larger the entropy increase. So the streamline 
passing through the nearly normal portion of the shock will have larger entropy 
increase than a streamline passing through the weaker portion of the shock. This 
change in entropy production in the nose region generates a layer of flow with 
entropy gradient, known as the entropy layer, which “influences” the development 
of the boundary layer in two different ways. First, it causes a continuous change 
of the flow properties at the outer edge of the boundary layer in the streamwise 
direction and second, it produces a velocity and pressure gradient at the outer edge. 
These external gradients invalidate the classical boundary layer approach. However 
Ferri [48] pointed out that these gradients are important when the vorticity of the 
inviscid flow is of the same order as that of the boundary layer. These conditions 
may exist when the Reynolds number is low and the Mach number is high. The 
thickness of the entropy layer is a function of the leading edge bluntness; even a 
small leading edge thickness generates an entropy layer which influences the local 
flow conditions for a large distance downstream. For some distance downstream of 
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the blunt nose the boundary layer grows inside the entropy layer which eventually 
is swallowed by the boundary layer in the case of flow past blunted cones. This 
distance is called the entropy swallowing distance. For conical flows, the entropy 
swallowing distance is defined [23] as “the location at the cone frustrum where fluid, 
which has gone through the strong portion of the bow shock, has been swallowed 
by the boundary layer.” This is by no means a precise definition as it depends 
upon the shock shape, boundary layer assumptions, and definition of the entropy 
layer thickness. The presence of the entropy layer has an important effect on the 
aerodynamic heating predication. 

For hypersonic flows over a blunted slender cone, the curvature of the shock 
produces the entropy layer as shown in Fig. 3.1. Once the shock attains conical 
shape, it is no longer curved and hence no more entropy gradient is produced. 
The thickness of the entropy layer decreases as it moves downstream on the body 
and eventually it is swallowed by the boundary layer. The entropy effects in the 
boundary layer vanish asymptotically further downstream. 

To calculate the entropy swallowing distance without calculating the flow 
field, certain assumptions regarding the shock shape and the boundary layer have 
to be made. Several investigators have attempted to approximate the shock shape 
for slightly blunted cones. Klaimson [49] has given one of the most accurate, yet 
simple, expression for the shock shape as 


y = 1.424cos0 c 



(3.1) 


For hypersonic speeds, the drag on the hemispherical nose can be calculated as 


C DT = 2 - cos 1 & c 


(3.2) 


A comparison of this shock shape with numerically calculated shape has been shown 
to give excellent agreement [49] within the region of interest. The value of y , where 
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the shock approaches the conical angle, is given by 


y c = cos 9 e 


0.984Cdt 


(sin 2 6 c ) 


0.426 


(3.3) 


The angle 6 C can be obtained from the relation 


Moo sin 6 e = 4 + 1.01 (Moo si 11 — 3.43) 


(3.4) 


With the known value of |7 C , the entropy swallowing distance, S e , can be approxi- 
mated as [50] 


1.5y* e M a 
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In Eq. (3.5), f(r/) is the transformed stream function defined as 

tP(S, V ) = {2§)$f(r,), /' = ^ 


(3.5) 


(3.6) 


where 


n = 
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I Jo Pe 


(2 S)i 
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In the present study, entropy swallowing distance is calculated by consid- 
ering the mass flow rate. The entropy layer thickness is calculated in the down- 
stream region, and the point where the entropy layer thickness becomes equal to the 
boundary layer thickness is termed as the entropy swallowing point. To estimate 
the thickness of the entropy layer, its edge is defined as the streamline that passes 
through the point A (Fig. 3.1) at which the shock attains the conical angle as given 
by Eq. (3.4). Most of the entropy changes occur before this point. The rate of mass 
flow in the entropy layer is given by 

rh = 2ir f puydy (3 7) 

J vt 

At point A, the mass flow rate is expressed as 

rh A = ny 2 e PooUoo 


(3.8) 
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Since the mass flow rate through the entropy layer is constant, y«(x) can be calcu- 
lated using Eqs. (3.7) and (3.8). At point A, y t is same as y c given by Eq. (3.3) 

The entropy swallowing distances calculated using present approach were 
compared with the results of Rotta [50] for various nose bluntness and Reynolds 
number. As will be shown later the present approach for calculating the entropy 
swallowing distance agrees qualitatively with the merging distance for skin friction, 
surface pressure, wall heat transfer and shock shape. The theoretical calculations of 
Rotta assume that the pressure on the conical portion of the body is constant and, 
therefore, the overexpansion and recompression on the shoulder are neglected. This 
causes the overestimation of the entropy swallowing distance. For small bluntness, 
the overexpansion can be neglected. Thus, Rotta’s results are valid only for very 
small values of nose bluntness. 

3.2 Governing Equations 

To calculate the flow over blunted slender bodies, two sets of equations 
are utilized. For the nose section, the use of full Navier-Stokes equations is neces- 
sary due to the presence of the subsonic region; but slightly downstream of the nose 
region on the body, the flow becomes predominantly supersonic and the Navier- 
Stokes equations are simplified for space marching. The Navier-Stokes equations 
without body force and external heat addition can be written in nondimensional, 
strong conservation form as [51] 

dU d(E - E v ) , d{F - F u ) l d{ G-G v ) = n (3g) 

aT + — ai +— Si az 

T 

U = {p,pu,pv,pw,pe t } 

E = {pu, pu 2 + p, puv,puw, (pe t + p)u} T 


where 
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and 


T 
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(3.10) 


p = (' y — l)p« (3-11) 

These equations are nondimensionalized with respect to the freestream 
values. The coefficient of viscosity p is obtained using the Sutherland’s law. For 
perfect gas calculations, Prandtl number is taken as 0.72. 

For the spherical nose region of the body, the flow is analyzed by the code 
SOFIA [52] . It solves complete Navier-Stokes equations using an explicit, finite 
volume, alternating two step scheme and uses local time stepping to accelerate the 
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convergence to steady state. It also incorporates time dependent grid adaptation to 
properly resolve the flow gradients without increasing the grid points. Once the flow 
field is obtained over the nose section, it is used as initial plane solution to initialize 
the parabolized Navier-Stokes equations. In the present investigation, a parabolized 
Navier-Stokes code developed by Vigneron et al. [53] and subsequently modified by 
Gnoffo [54] is used for analyzing the downstream flow over the bodies. Parabolized 
Navier-Stokes (PNS) equations are subset of the complete Navier-Stokes equations 
and can be used to predict complex three-dimensional, steady, supersonic, viscous 
flow fields. These equations can be solved using a space marching technique as 
opposed to the time marching, which is usually employed for the Navier-Stokes 
equations. They are valid in both viscous and inviscid regions and thus unlike 
boundary layer equations no special effort is needed for viscous/inviscid interactions. 
They can be obtained by 

(i) Neglecting unsteady terms 

(ii) Neglecting streamwise viscous diffusion terms 

(iii) Modifying the streamwise flux to permit the stable time like marching of the 
equations from the initial data plane. . 

For these equations to be valid, the flow in the inviscid region should be supersonic 
and stre am wise velocity component should be positive. Thus, these equations fail 
in the presence of streamwise separation, although they are valid in the region of 
crossflow separation. The streamwise pressure gradient needs special treatment as 
will be explained later. 

The following two independent variable transformations are applied to al- 
low for conical effects and stretching between the body and shock. 

c = - (3.12) 

X 

f = f(6,c) (3.13) 


a = x 


6 = - 


| = a 


rj = r?(a,6,c) 
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Upon combining Eqs. (3.9), (3.12) and (3.13), neglecting unsteady terms and vis- 
cous derivatives in £ direction, the PNS equations are expressed as [53] 

(3.14) 


dE x dF\ dGi _ q 

d£ dr] d$ 


where 


Ei = 


a 2 E 


T djrht) 

d(b,c) 

The presence of streamwise pressure gradient term allows for upstream 
propagation of information in the subsonic regions of the flow; thus, the equations 
are not well posed for space marching method. Various remedies have been proposed 
to overcome this problem, the most successful one was proposed by Vigneron et 
al. [53]. In this approach, a fraction of the pressure gradient term is retained in the 
subsonic region and the remainder is either omitted or is evaluated explicitly using 
the backward differencing technique. Thus the term E\ is split as 

E x = E* + P ( 3 - 15 ) 

T 

E* = {/>u,pu 2 + u}p,puv t puw,[pc t + P)u} 

P = {0,(1 -w)p, 0,0, 0} T 

This treatment of the pressure gradient allows for stable marching scheme without 
encountering departure solutions. An eigen value analysis given in [53] indicates 
that for stability 


„ < (3-w) 

1 + (t - 1)MJ 

where w is the factor of safety; in the present calculations its value is taken as 0.9. 
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3.3 Method of Solution 


The PNS equations are a mixed set of hyperbolic-parabolic equa- 
tions in the streamwise direction. Equation(3.14) is solved with a finite difference 
method using the Beam and Warming algorithm [55]. It is an implicit, non-iterative 
algorithm with approximate factorization in Delta form. The Vigneron technique 
is used to suppress the numerical instabilities due to streamwise pressure gradients. 
Details of n um erical procedure are available in [53]. For the sake of completeness 

only a brief description is given here. 

Using the Euler implicit scheme, Eq.(3.14) can be written in the Delta 


form as 


' aE r + A4(^)id£r‘[|£ + Ae£(f£)]A‘p. 


aUi -'d('du l '"av l ’ Wi 


dri'dU,’ 


(3.17) 


where 

Ih = a 2 U/J, A*U! = U\ +1 - U[ and A e P = A‘ -1 P 

Since the vectors E\, F u and Gi are homogeneous functions of degree one in U, 
the conservative form of the governing equation is maintained. The viscosity n is 
assumed independent of vector U and is a function of the position only. The cross 
derivative ter ms are neglected in the calculation of Jacobian. Equation(3.17) is 


solved in the following four steps. 


i|f +A 40 Aipr=MS(3 - 17) 

A*-$f>ATO 

[2S + A^— (— )]A‘f7i = AZTT 
[ dU x Z dri K dlh n 

U{ +1 = U{ + A U[ 


(3.18) 


This procedure avoids the computation of inverse matrix and introduces two one- 
dimensional operators which are solved using the Thomas algorithm. In order 
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to d am p the spurious oscillations, second order implicit and fourth order explicit 
damping terms are added. After adding the dissipation terms the truncation error 
of the algorithm is consistent with the first order Euler scheme. 

Due to symmetry of the flow field, only one half of the flow field is com- 
puted. The generalized coordinates rj and £ are defined such that the computational 
plane has a square shape of side unity with uniform spacing in both directions. The 
outer boundary corresponds to the shock and the inner to the body. The grid 
points are clustered near the body to properly resolve the gradients. A typical grid 
is shown in Fig. 3.3. 

At the body surface, no slip, zero pressure gradient and constant wall 
temperature conditions are imposed. At the outer boundary, shock fitting approach 
is used and the flow variables behind the shock are calculated by the Rankine- 
Hugoniot relations. At the plane of symmetry, reflecting conditions are used. As 
mentioned earlier, the initial plane solution to initiate the marching procedure is 
provided by the Navier-Stokes equations. 
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Fig. 3.3 Typical grid for finite difference calculations. 
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3.4 Equilibrium Chemistry 

To include the real gas effects which might be present at hypersonic 
speeds, calculations are also made with equilibrium air chemistry. The equilibrium 
chemistry package developed in [56] is used in the analysis. In this model, complete 
thermodynamic equilibrium is assumed, in which the gas properties and concentra- 
tion of atoms and molecules are identical with their equilibrium values appropriate 
to the local conditions. The rates of reaction are assumed very fast in comparison 
with the rate of diffusion. Here, the free energy minimization technique with the 
method of steepest descent is utilized. The mole fractions of the constituent species 
are calculated as a function of temperature and pressure. The assumption of equi- 
librium is made on the basis of physics of the problem. In the blunt body flows 
near the stagnation point, where the shock is nearly normal, the velocities are very 
low and temperature is very high. Consequently the flow attains equilibrium values 
prior to leaving the nose region. The particles crossing the bow shock farther from 
the axis experience higher velocities and lower temperature due to weaker shock, 
thus reducing the degree of non-equilibrium. On the after body past the sonic line, 
the flow experiences expansion thus increasing the velocity while reducing pressure 
and density, both of these effects try to quench the chemistry and produce frozen 
flows. Thus the afterbody flows can be calculated by using an effective value of 7 
determined from the mas s concentration. Similar approaches have been taken pre- 
viously by several authors [57,58]. The equilibrium chemistry option can easily be 
incorporated into a perfect gas code with equivalent 7 approach where 7 I s now de- 
fined as the ratio of enthalpy to internal energy. Further details on implementation 
of the equilibrium chemistry model are given in Appendix A. 
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3.5 Validation 


The PNS code was validated by comparing the present results with 
those calculated using the laminar theory with reference temperature concept [59]. 
For the same Reynolds number, the skin friction on a cone and a flat plate is related 


by 


C ft, cone ~ VSC 


fe, plate 


Consequently, the heat transfer coefficient is expressed as 


(3.19) 


^ fe,cone 

'he, cone — “ _ 2 

2 Prt 


The value of Cf e , p i at e can be calculated from 

0.664\/C* 


'fe,plate 


(3.20) 


(3.21) 


where Re xt is the Reynolds number based upon the distance along the cone surface 
and the subscript e refers to the edge conditions. The surface flow along the cone is 
assumed to be constant, then p e , U e , and T e can be calculated by inviscid supersonic 
cone theory [60]. The quantity C* in Eq.( 3.21) is the reference viscosity and is given 
by the Sutherland’s law as 


C* 


/T *\ 06 T e + K 
V7W T* + K 


where 

K = 200° iE 


and 

= 0.5 + 0.039M* + 0.5^ 

T t ±e 

By using Eqs.(3.19)-(3.21), the theoretical values of Cf t%con e and Che, cone 
calculated. 

These results are compared with the numerically calculated values for 
sharp cone at A/qo = 10 and are discussed in Chap. 5. 



Chapter 4 


EFFECTS OF SHOCK-SHOCK 
INTERACTIONS 


In this chapter, the physical models used to calculate the two-and three- 
dimensional shock-shock interactions in the scramjet inlet are discussed. The vari- 
ous types of interaction patterns as described by Edney [31] are reviewed and how 
they occur in the scramjet configuration is discussed. The governing equations and 
method of solution used are described. Depending upon the Mach number and the 
angle of attack, the forebody shock can interact with the inlet flow; thus producing 
a very complex flow field. The primary factors affecting the nature of the interaction 
are freestream conditions, strength of the impinging shock, geometry and impinge- 
ment location. Shock-shock interaction can also occur on space shuttle, missile 
launcher and other high-speed aircrafts. As mentioned earlier, this study concen- 
trates on interactions in scramjet inlet but the results obtained are very general in 
nature and can be applied to other applications as well. All interactions considered 
here are of the free interaction type i.e., the region upstream of the interaction 
zone is independent of the source of impinging shock. Two types of interaction can 
occur in the scramjet inlet; interaction of the forebody and the cowl shock, and the 
interaction of the forebody and the inlet sidewall shock. The interaction of cowl lip 
shock with the forebody shock is modeled by a cylindrically blunt wedge with an 
oblique shock impinging on it in such a way that the line of interaction is parallel 
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to the cylinder axis. If the end effects are neglected, the flow in each plane normal 
to the interaction line should be identical. Therefore, a two-dimensional model as 
is shown in Fig. 4.1 is used. The interaction of inlet sidewall shock and the fore- 
body shock is more complicated because the inlet sidewalls are swept back and the 
forebody shock is normal to the axis. Thus the flow is fully three-dimensional in 
nature. The interaction point of the forebody and inlet sidewalls shocks is far away 
from the inlet/body junction. Therefore, the effect of the body on the shock-shock 
interaction is neglected and the inlet sidewalls is modeled by an infinite swept blunt 
wedge with impinging shock normal to its axis as shown in Fig. 4.2. 

4.1 Classification 

Edney [31] in 1968 categorized the shock-shock interactions into six basic 
patterns (Types I to VI, Figs. 4.3 - 4.8). The primary factors in determining the 
type of interaction are strength of the impinging shock, location of the impingement 
and the freestream Mach number. All six types of interactions can occur when an 
oblique shock interacts with the bow shock of the leading edge. This classification is 
a general one and sometimes it is not clear to which category an interaction belongs. 
Types I, II, and V interactions result in shock-boundary layer interactions, Type 
III results in attaching shear layer interaction, Type IV produces a supersonic jet, 
surrounded by subsonic flow and the Type VI interaction results in an expansion 
fan-boundary layer interaction. The shock-shock interaction causes the stagnation 
point to move and produce a localized zone of high pressure and heat transfer. 
Since for a blunt body the peak heating and pressure occur at the stagnation point, 
an amplification factor is defined as ratio of the peak value of pressure (or heat 
transfer) divided by its value at the stagnation point for the unimpinged case. 

The Type I interaction occurs when two weak shocks of opposite family 
interact. When an oblique shock intersects the leading edge shock well above the 
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Fig. 4.2 Physical model for forebody and inlet shock interaction. 
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upper sonic line, the interaction occurs as is shown in Fig. 4.3. The interaction 
produces a shear layer and a transmitted shock; the shock strikes the body, thus 
causing shock-boundary layer interaction. If the transmitted shock is strong, it can 
cause boundary layer to become turbulent or even separated. Holden [61] correlated 
the peak pressure and heat transfer with the relation q/q 0 = (p/p<,)° 85 over the Mach 
number range 2.4 - 13. 

A Type II interaction occurs when two shocks of opposite family intersect. 
If the impinging shock intersects the leading edge shock just above the upper sonic 
line, this type of interaction is produced. This interaction produces two shear layers 
and a transmitted shock as shown in Fig. 4.4. The transmitted shock strikes the 
body, causing pressure and heating amplification. Depending upon the body shape, 
the shear layer can also strike the body, thus resulting in additional amplification 
of heating and pressure. In general, the upper shear layer(from pt. A) is stronger 
than the lower shear layer ( from pt. B). 

A Type III interaction occurs when a weak shock wave interacts with the 
strong shock wave. When the oblique shock strikes the bow shock inside the sonic 
region, this type of interaction occurs. The flow is supersonic above the shear layer 
and subsonic below it. A reflected shock exists between the body and the leading 
edge shock which turns the flow parallel to the body. The pressure and heating 
rates are amplified due to striking sheax layer and reflected shock, times. 

A Type IV interaction occurs when the impinging shock strikes the lead- 
ing edge shock near the stagnation zone, where the shock is nearly normal. This 
produces a supersonic jet bounded by two shear layers and submerged in subsonic 
flow. Near the body, it produces jet bow shock and stagnation zone when it strikes 
the body. This produces a very complex flow field with presence of shocks, shear 
layers and the jet. Type IV interaction is the most severe case and produces the 
largest amplification of heating and pressure. 


Fig. 4.3 Type I 
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-shock interaction. 







Fig. 4.6 Type IV shock-shock interaction. 
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Fig. 4.7 Type V shock-shock interaction. 
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Fig. 4.8 Type VI shock-shock interaction. 
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When the oblique shock intersects the bow shock just below the lower 
sonic line (Fig.4.7), Type V interaction is produced. In this case, the shocks are 
of the the same family. The interaction produces a thin supersonic jet and a shear 
layer. It is very similar to the Type II interaction. 

A Type VI interaction occurs, when the two weak shocks of the same 
family interact; these shocks are weaker than Type V interaction. Such a situation 
exists when a weak impinging shock intersects the bow shock well below the lower 
sonic line. An expansion fan is produced from the intersection point, which causes 
decrease in pressure and heating rates. 


4.2 Governing Equations 


The Navier-Stokes equations in cartesian coordinates describing three 
dimensional, unsteady compressible flow in the absence of external body force can 
be written in differential form as [51] 


dQ dtfjJQ d(F - F„) a(G-G.) _ 0 

dt dx dy dz 
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The molecular viscosity is calculated from the Sutherland’s law and the Stoke’s 
hypothesis for bulk viscosity, A + ^ = 0 is invoked. The perfect gas equation of 
state is used to define the pressure p via the internal energy as 


P 


(7-1) 


t — p(u 2 4- v 2 + tu 2 )/2 


(4.2) 


4.2.1 Nondimensionalization 

Equation(4.1) is nondimensionalized as follows 
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In the above equations tilde denotes the dimensional quantity and the subscript oo 
denotes the freestream quantity. 



4.2.2 Coordinate Transformation 


The governing equations are transformed to the body-fitted coordinate 
system (£, r],$) using a steady transformation of the type 


Z = Z{x,y,z) 

r) = n{x,y,z) 

$ = ${x,y,z) 


The Jacobian of the transformation is denoted by J and is given by 

j _ d(g,T7,f) 

d(x,y,z) 


x € y„ 2 f + x { y ( z n + x r ,y i z i - x ( y ( z v - x n y ( z t - x ( y n z ( 


(4.3) 


Using the chain rule of partial differentiation, the transformation metrices are de- 
fined as follows [51] 


tx/J = y, ,z ( - z n y { 
Zy/ J = z r)X( X r jZ( 

Zu/J =.Xr,y t - y n x i 
r) x /J = z ( y ( - y ( z { 
Vy/J = x iZ( - Z(X( 
Vz/J = Vi x ( - x (y ( 

Sz/J = y^Zr, - z&r, 

-- Z^Xr, - X(2„ 

$z/J = n - yt x n 


After applying the generalized coordinate transformation to Eq.(4.l), it can be 
expressed in the strong conservation- law form as 
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where the new variables are defined as 
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The inviscid and viscous flux vectors in generalized coordinate system can now be 


defined as 
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where U, V, and W are contravarient velocities and are given as, 


U = £*u + £„v + 
V = r? z u + r7 v v + 
IV = 
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and 

K = u i r *.*y - 9*. ( 4 - 5 ) 

For two-dimensional calculations, E and E v were set equal to zero, thus reducing 
the governing equations to two-dimensions with f normal to the surface and along 
it. For three-dimensional case, viscous terms were dropped in £ and rj directions 
due to coarse grids in axial and circumferential directions. 

4.2.3 Finite Volume Formulation 

Equation ( 4.4) is solved using a finite volume method. In this method, the 
integral formulation of the conservation laws is discretized directly in the physical 
space. The direct discretization ensures that mass, momentum and energy are 
conserved at discrete levels. It remains valid in the presence of discontinuities in 
the flow field such as shocks and contact surfaces. Eq.( 4.4) is expressed in the 
integral form as 

§iIIL QdV+ Ils phdS=0 (4 - 6) 

where F = (E - E v )i + (F - F v )j + (G - G v )k, and h = n x i + n v j + n z k is the 
unit normal vector pointing outward from the surface S bounding the volume V. A 
semi-discrete finite volume representation of Eq.(4.4) leads to 

+ (£- -(e- 

+ (F — Fv)ij+l t k — {F ~ 

+(<j - G.hjMi - < 6 - = 0 ( 4 - 7 > 

It can be shown that Eq.(4.6) can be written in a form identical to Eq.(4.7) by 
writing the surface integrals as the sum of the contribution from six faces of the 
hexagonal cell. In the computational domain(£,77, f), 
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A£ = e i+ i-6_i = i 

Similarly, Ar? = Af = 1 

The fluxes are evaluated at the cell boundaries defined by the grid points 
as is shown in Fig. 4.9. So far the method of solution has been identical to the 
finite difference method. In most of the finite difference methods the transformation 
matrices and Jacobian is calculated by using the two point central difference formu- 
las. This leads to inconsistency in the volume and surface normal calculations such 
that the geometric conservation laws are not satisfied and hence the scheme fails to 
capture the freestream when evaluated on arbitrary meshes. In the present method, 
geometrical interpretation of the Jacobian and transformation terms is made. The 
Jacobian is calculated as inverse of the cell volume, the vector A k/ J is the directed 
area of the cell surface to a k = constant coordinate direction and \Ak\/J is the 
area of the cell interface. Here k represents £,»7, and f, respectively. 

The volume of a cell is defined by its eight vertex points and is calculated 
as the sum of six pentahedrons. Each pentahedron is formed by one of the six cell 
faces and the average point in the cell volume. The directed area of the cell interface 
is calculated as one half of the vector cross product of two diagonal vectors joining 
the four vertex points. In general, all four vertex points do not lie in the same 
plane. So the directed area so formed corresponds to the planar surface obtained 
by passing a least square plane through the four points and then projecting them 
into the plane. This approach satisfies the geometrical conservation laws and makes 
the numerical scheme compatible with the finite volume formulation. 

4.3 Method of Solution 

Over the past few years, upwind space discretization methods have be- 
come a popular alternative to the central space discretization methods. The basic 
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Fig. 4.9 Finite volume formulation. 
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idea behind the development of upwind scheme is that they try to mimic the direc- 
tional dependence by applying asymmetric space discretization with a bias toward 
the direction from which the information is propagating. The upwind differencing 
reflects the predominant nature of the Navier-Stokes equations in the limit as Re 

— oo (hyperbolic) and Re -*• 0 (parabolic). This improvement in physical treat- 
ment comes with increased computational work per iteration but the total number 
of iterations is reduced thus offsetting the increased computation time per iteration. 
These methods are naturally dissipative and therefore no artificial viscosities with 

problem dependent coefficients are needed. 

The Navier-Stokes equations consist of inviscid (pressure and convective) 
terms and viscous (diffusive) terms. In the upwind methods, the inviscid fluxes 
are split into positive and negative components based on eigen values, and then 
spatially differenced by backward or forward differencing. The differencing could 
be first order, second order fully upwind or third order upwind biased. First order 
upwind differencing have a large amount of numerical dissipation which makes it 
unsuitable for meaningful calculations. The second order upwind schemes have 
second order truncation error; the leading term of this truncation error is dispersive 
and the leading numerical dissipative term is actually a third order, fourth derivative 
term. Even for third order differencing the scheme is actually second order because 
of the second order treatment of diffusion terms. Hence, the numerical dissipation 
of the second order upwind scheme is of the same order as the central difference 
scheme. The central difference schemes have been shown to be superior than upwind 
schemes in the viscous zone [62]. 

Various upwind algorithms are available in the literature. These are A 
scheme [63], split coefficient method (SCM) [64], flux difference methods [65] 

- [68] and flux vector differencing [69] and [70]. For the physical problem of shock- 
shock interaction, various schemes were tried. The first scheme tried was explicit 
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MacCormack [71] with central differencing. The solution had excessive oscillation 
at high Mach numbers and required proper tuning of damping coefficient, also due 
to explicit nature of the scheme, the time step was very restrictive for the grid 
employed. The flux differencing algorithm developed by Roe [68] gave physically 
incorrect solution even for unimpinged blunt body flows. The bow shock had a 
dimple near the stagnation point. For some cases, the shock did not even stablize 
but continued to move upstream until stopped by the inflow boundary. The peak 
pressure value was not at the stagnation point but rather away from it. The same 
type of phenomenon has been observed for the Roe’s scheme by Perry et al. [45]. 
The critical Mach number around which this phenomenon starts is about 6.0 [72]. 
This phenomenon is not fully understood yet, and therefore no rigorous correction 
to Roe’s scheme to eliminate this problem is available. However, some progress has 
been made for simple flows by adding eigen value smoothing [45]. van Leer’s [70] 
flux vector splitting method gave the best results for this problem. Consequently, 
this method has been used for two-and three-dimensional shock-shock interactions. 

4.3.1 Flux Vector Splitting 

Steger and Warming in 1981 [69], realizing the homogenity property of the 
Euler equation, proposed a novel way of splitting the inviscid fluxes into positive and 
negative subvectors based on the eigen values. The positive subvector is associated 
with the non-negative eigen values and negative subvector is associated with the 
non-positive eigen values. The split fluxes are then spatially discretized using one 
sided upwind differencing. The splitting is not unique and various combinations 
are possible as long as the sum of two subvectors is the same as the original vector. 
The details of this splitting and its implementation into Euler equation is available 
in [69]. The main drawback of this type of splitting is that forward and backward 
fluxes although continuous, axe not differentiable when an eigen value changes sign. 
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Such phenomena occurs at stagnation and sonic points. This causes small glitches in 
the solution near these points, van Leer in 1982 [70] proposed a similar splitting but 
with continually differentiable fluxes by introducing additional requirements that 
the Jacobian must be a continuous function of the Mach number. As mentioned 
earlier, upwind differencing is applied to inviscid fluxes only, the viscous fluxes are 

centrally differenced. 

Splitting the fluxes in £ direction, 

E = E + + E~ ( 4 - 8 ) 

where 

E + = E , E~ = 0 for > +1 

E~ = E , E + = 0 for < -1 

Here is the contravariant Mach number along £ and is given as 

M ( = v/c , V = U/\At\ 

For subsonic flows. |Af*| < 1, the split fluxes are obtained as follows 
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The conserved state variables on the cell interface are obtained as 

Q ~ + 1 = Qi + (1 — K f) V* +(1 + 

Q, + + , = Qi+i - ^«[(i + *e) V £ +(1 - /c { )A«]g <+1 


(4.10) 
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where 

A (Qi = Q<+i — Qi 
ViQi = Qi~ Qi - 1 

The switch <j> is zero for the first order differencing and unity for higher order differ- 
encing. The parameter k determines the order of differencing; = — 1 corresponds 
to the fully upwind second order, = 0 corresponds to second order upwind bi- 
ased, = +1 to second order central differencing, and k ? = +1/3 to the third 
order upwind biased. Having obtained the value of Q at cell interfaces, the split 
flux differences are implemented as follows 

6 ( E, = 6 ( E~ + 6 ( E; ! 4 - u ) 

= £♦(«-+ i - i+e-li-j + E- w*) i+4 - E- (<3 + ),. t 

The term j& + (Q“), + ^ denotes the forward flux evaluated at the cell interface t + 
I. it should be noted that MUSCLE ( monotonic upstream-centered scheme for 
conservation laws) approach is used. The conservation variable Q is obtained first 
at the cell interface and then the split fluxes E ± are obtained. This is in contrast 
to obtaining E ± first at the cell center and then interpolating it to get the cell 
interface value. The former approach has been shown superior than the latter as 
the split fluxes are less differentiable than conserved variables when transitioning 
through sonic and stagnation points. Similarly, the fluxes can be split in tj and f 
directions also. 

4.3.2 Algorithm 

For the present calculations, the algorithm developed by Thomas et 
al. [73] - [75] has been used. The algorithm is outlined for three-dimensional flows 
and can be specialized for two-dimensional flow as necessary. As noted previously, 



viscous terms have been dropped in £ and rj direction; thus Eq.(4.4) can be written 


as 
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Applying the Euler implicit time integration scheme and linearizing it, one obtains 
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Equation(4.12) represents a large banded block 5x5 matrix equation, the 
band width depends upon the grid size and choice of spatial discretization. Using 
second order accurate spatial discretization, it can be expanded as 


(jZi + C( + C " + C ' )AQ " + 

-tA t AQ”_, + B(AQ"_ 1 + D ( AQ?+, + H C AQ%, 

+A tJ AQ'j_2 + BqAQ"-! + D r} AQ r j +l + H n AQj +2 

+A ( A<5*_2 + B{AQk_i + DfAQI+i + H { AQ k+2 — R{Q ) (4-14) 


where A,B,C,D and H are 5x5 block matrices associated with implicit spatial differ- 
encing in the £,q, and f directions. The solution of Eq. (4. 14) involves inversion of 
a block pentadiagonal matrix equation. The computational effort is reduced if the 
implicit-spatial discretization is taken as first order accurate. This does not affect 
the steady state accuracy which is determined by the spatial differencing of R(Q) • 
Then, A € , A ( ,H^H„, and are equal to zero and 
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c s = (415) 

_ , d b-\ 

D *~ ^ dQ ' ,+ 1/2 

Similar expressions can be developed in rj and f directions also. 

Using Eq.(4.16), Eq(4.14) can be written as 

( + C( + C, + C,)AQ" + B ( AQ?_, + D t AQ; + , 

+ R.AQ'V, + B, A<?;_, + D s AQ" m = R(Q") (4-16) 

Relaxation is implemented in the streamwise direction and the approximate factor- 
ization is implemented in the cross flow planes. 

For first order implicit space discretization, the algorithm is stable regard- 
less of the sweep direction. With higher order spatial differencing, alternate forward 
and backward sweeps are required. For a forward relaxation sweep, AQ 1 n _ 1 is known 
and AQ" +1 = 0, then Eq.(4.16) becomes 

(yAj + <7, + C„ + C,)Ag“ + a.Ag;., + D,A<?7 +1 
+B-aq;_, + £>,a <?; +1 = «?«.«?«) < 417 ) 


Equation(4.17) is written in a compact form as 
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Equation (4.18) is approximately factored in rj and f directions as 

= r(q?-2, q?, q? + i , q: + 2) 


(4.19) 
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The solution of Eq.(4.19) then involves the solution of two block tridiagonal equa- 
tions 

r)F 

\M + 6,^ = R(Q?h Q?, Q? + QU) 

gn+l _ Q» + (4.20) 


For the cowl afterbody calculations the algorithm has been modified so that 
the downstream influence of the subsonic part of the boundary layer in an otherwise 


supersonic flow is suppressed by restricting the streamwise pressure gradient. A 
parabolized solution is then obtained by marching downstream and iterating locally 
in each plane until convergence. The Vigneron [53] technique as described previously 


in Chap. 3 is used to modify the algorithm. The details of modified algorithm are 
available in [74]. The forward flux vector E is now modified as 
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where u is the Vigneron factor and was defined in Eq. (3.16). The additional terms 
due to the explicit lagging of the streamwise terms are neglected, then the procedure 
is not subjected to the departure solution for arbitrary streamwise mesh spacing. 
The modified flux and flux Jacobian are split as 


E + {Q )i± 1/2 — &{Q)i± 1/2 (4-22) 

d£ + (Q-) = dE{Q) 
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and 
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(4.23) 
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The solution from the (i-1) station is used as initial guess at ith station and is 
iterated locally to obtain the converged solution and the process is repeated at the 
next (i+1) station. The resulting algorithm is fully second order accurate. The 
value of k is taken as -1 in order to have fully upwind algorithm. 

In order to eliminate the overshoots and undershoots near the shocks, 
Van Albada type continuous differentiable flux limiter has been used in the present 
calculations. It reduces the scheme to fully one sided. The limiters modify the 
upwind biased interpolation and are implemented by modifying Eq.(4.10) as follows 


Q i+ i -Qi + 4 


3, + + * - Gi+i 4 |_ 


(1 - SAC*) V{ +(1 + S/C()A{]<2,- 

(1 + s/c £ ) V« +(1 - «« £ )A £ ]q, +1 


(4.24) 


where 

^ _ 2A t V( + { 

S “ A| + V| + « 

and e =10 -6 is a small number to prevent the division by zero in the region of zero 
gradients. 


4.4 Grid Generation 


The first step in any numerical solution is the discretization of the gov- 
erning equations from a continuous domain to a set of discrete points (called grid 
points). The choice of grid points is not an arbitrary one but is governed by physics 
of the problem. Grid generation is a procedure for orderly distribution of the sam- 
pling stations over a physical field in such a way that all physical phenomena in the 
entire region of interest may be represented with sufficient accuracy by this finite 
number of sampling stations. Since the computer memory and speed limit the num- 
ber of grid points, it is very crucial to make best use of the available resources. The 
grid points should conform to the boundaries and be concentrated in the regions 
of high gradients such as shocks and boundary layers in order to properly resolve 
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the flow field. They must also be oriented in such a way that grid is as close to 
orthogonal as possible. 

The grids can be generated by either numerically solving the partial differ- 
ential equations or algebraic equations. In the present study, an algebraic method 
is used. In this method, the grid is generated by algebraic interpolation from the 
boundaries. Here an explicit functional relationship between the computational and 
physical domain is known. It is much faster than the differential equation approach 
and offers an easy control of grid spacing and distribution. For the present calcu- 
lations, various types of grids were employed. For the first set of calculations, an 
outer boundary consisting of a circular arc was defined. Since the shock capturing 
algorithm has been employed, the outer boundary was chosen in such a way that 
there were sufficient number of points between the shock and outer boundary. The 
grid points were clustered normal to the body using the exponential stretching 


to resolve the boundary layer. Since the heat transfer and skin friction are strong 
functions of grid spacing, it is very important to properly resolve the gradients. For 
the impingement case the shock moves closer to the body on one side and away from 
the body on the other. Thus a lot of grid points are wasted. In order to alleviate 
this problem, a semi adaptive grid was used. Figure 4.10 shows this grid, in this 
case the outer boundary is made to conform with the shape of distorted bow shock, 
which was obtained by coarse mesh calculations where no attempt was made to 
resolve the gradients. This type of grid, reduced the wastage of grid points but it 
requires a new grid every time either freestream conditions are changed or the type 
of interaction changes. The solution obtained by using this type of grid was much 
better resolved near the shock region but in the case of impingement where we have 
a supersonic jet enclosed within subsonic flow and extremely high gradients, the 
resolution was still less than desired, as will be shown in Chap. 5. 
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In order to improve upon the solution, adaptive grid system has been 
used, which adapts the grids during the course of the solution in order to follow 
the developing gradients in the physical solution. The grid points move as the 
solution develops, concentrating the points where they are needed the most. The 
total number of points were kept constant. The final adapted grid for Mach 8.03 
case is shown in Fig. 4.11. For this, the method developed by Abolhassani et al. [76] 
has been used and is explained in Appendix B. It is a very general method with 
capability to adapt the grids with various variables such as pressure, Mach number, 
shear stress etc. and is based upon the variational approach. Since the solution 
varies predominantly in the normal direction, the grids were adapted in only one 
direction. 

The grid for three-dimensional calculations is shown in Figure 4.12. The 
grid points were clustered near the body in normal direction as well as near the 
shock impingement location to resolve the gradients. As in the previous case, the 
outer boundary was moved far enough so as to avoid any interference between the 
shock and the outflow boundary. 
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Fig. 4.10 Unadapted grid for two-dimensional calculations 



Fig. 4.11 Adapted grid for two-dimensional calculations 



Fig. 4.12 Grid for three-dimensional calculations. 




Chapter 5 

RESULTS AND DISCUSSION 


The numerical schemes described in Chaps. 3 and 4 are applied to obtain 
steady state results for the forebody and the shock-shock interaction on the cowl 
plate and inlet sidewalls. All calculations were performed under the assumption of 
a laminar flow. For the forebody, the flow is calculated by the parabolized Navier- 
Stokes equations with the starting solution for the nose section provided by the 
full Navier-Stokes equations. For two-dimensional shock-shock interaction, the full 
Navier-Stokes equations are used and for three-dimensional interaction, the thin- 
layer Navier-Stokes equations are used. Results for the nose bluntness effects on the 
forebody are discussed first followed by the results for the shock-shock interaction 
on the cowl plate and inlet sidewalls. 

5.1 Parametric Study of Nose Bluntness Effects 

The parametric studies were conducted to determine the extent of down- 
stream effects of nose bluntness for various freestream conditions given in Table 5.1 
for three different values of nose bluntness. The blunted cone with the smallest 
nose radius is considered as a sharp cone and the effects of leading edge bluntness 
are compared with respect to this cone. When the surface quantities fall within 
five percent of the sharp cone values, the effect of blunt ness is assumed to have 
vanished. The same criterion is used for the case of ogive also. 
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First, the results for a 5° sharp cone at Mach 10 are compared in Figs. 
5.1 - 5.4 with the results of inviscid and laminar theory to validate the numerical 
procedure as well as to justify the assumption of smallest nose radius cone as a 
sharp cone. Next the results for the wall quantities are presented for various Mach 
numbers and cone angles in Figs. 5.5 - 5.12. The effects of Reynolds number on 
the wall quantities are presented in Figs. 5.13 - 5.16. Similar results for ogives are 
presented in Figs. 5.17 - 5.23 for wall and flow field quantities. 

The numerical results are validated by comparing the surface pressure and 
shock standoff distance with the inviscid theory in Figs. 5.1 and 5.2. The pressure 
has been nondimensionalized with respect to the freestream pressure and the shock 
standoff distance with respect to the nose radius. The x distance is measured 
from the virtual tip of the sharp cone. The results are plotted from the tangency 
point onward. The skin friction and the heat transfer coefficient are compared 
with the laminar theory in Figs. 5.3 and 5.4. These figures show an excellent 
comparison. The theoretical results are valid for sharp nose cone only, while the 
numerical calculations were performed for smallest nose radius cone. Except very 
close to the nose, the results agree with the sharp cone results. So the assumption 
of smallest nose radius as sharp cone is valid slightly downstream of the nose. 

Figures 5.5 - 5.8 show the variation of surface quantities with the 
axial distance for a 5° cone with various degree of bluntness and for freestream 
Mach numbers of 10 and 20. The surface pressure distributions are presented in 
Fig. 5.5. In the stagnation region, the pressure is very high but as the flow moves 
downstream, the shock strength decreases thereby decreasing the pressure. The 
flow overexpands near the shoulder; the extent of overexpansion is a function of 
the bluntness and freestream conditions. The overexpanded flow then recompresses 
back to the sharp cone value. The distance it takes to reach the sharp cone value 
is also a function of the freestream conditions and nose bluntness. It should be 
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Fig. 5.4 Comparison of heat transfer coefficient with the laminar theory. 



Table 5.1 Flow conditions for the forebody calculations 


Quantity 

Moo = 10 

.* 

II 

to 

o 

Poo,N/m 2 

404.7 

171.0 

Too , k 

243.4 

261.3 

Tyjall > ^ 

1000 

1000 


Pr 


0.72 


0.72 
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noted that the pressure on the cone is influenced by the nose over a large portion of 
the afterbody and is lower than the conical pressure. For all values of leading edge 
bluntness considered, the surface pressures have reached the same value at about 
160 nose radii downstream. 

Figure 5.6 shows the variation of the shock standoff distance as a function 
of the axial distance. The shock standoff distances are measured from the body 
surface. Here the behavior of the shock shape is seen to be qualitatively similar 
to those observed by previous investigators [20]. The shock standoff distance is 
affected considerably by the leading edge bluntness. For the sharp cone, it is a 
straight line, while for blunted cones it is curved with inflection point about 50 
nose radii downstream. For sufficiently large values of x, the shock shape becomes 
independent of the nose bluntness. At Mach 20, it takes less distance to attain 
conical value than at Mach 10. 

The variation of the wall heat flux with axial distance is illustrated in Fig. 
5.7 for various freestream conditions and bluntnesses. The wall heat flux has been 
nondimensionalized by pooU^- Figure 5.7 shows that the heating rates decrease 
over the body with increasing nose bluntness at a given Mach number. This is 
because heat transfer in the nose region is inversely proportional to the square root 
of the nose radius and decreases over the afterbody due to lower local pressure. The 
results again show that the sharp cone value is attained in a shorter distance for 
Mach 20 than for Mach 10. Comparing the distance it takes to attain the conical 
value for pressure and heat flux, it is noted that it takes almost the same distance for 
both. Since (to a reasonable accuracy) the pressure can be predicated analytically 
by using various correlations, it appears that the merging distance for heat transfer 
can be estimated without calculating the detailed flow field. 

Variation of the skin friction coefficient with the axial distance is shown in 
Fig. 5.8. The skin friction is zero at the stagnation point, it increases around the 
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Fig. 5.6 Variation of shock standoff distance with the axial distance for 6 C = 5° 
various Mach numbers. 
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nose, reaches a maximum and decreases further downstream. It follows the same 
general trend as the wall heat flux over the afterbody. Figure 5.8 shows that the 
skin friction decreases over the body with increasing nose bluntness at a given Mach 
number. This is because the entropy layer increases the thickness of the boundary 
layer, thereby reducing the gradients of velocity profile near the body. The skin 
friction is found to be most sensitive to the nose bluntness as it takes maximum 
distance to reach the sharp cone value among all the flow properties. 

A comparison of the entropy swallowing distance with the results of Rotta 
[50] for Mach 10 freestream conditions (Table 5.1) is given in Table 5.2 for various 
nose bluntness and Reynolds numbers. The table shows a large difference between 
the present results and the empirical results of Rotta. This difference increases 
with the increasing bluntness. The calculated entropy swallowing distance agrees 
with the merging distance for skin friction, surface pressure, wall heat transfer 
and shock shape. As mentioned earlier, the discrepancy between the two is due 
to the assumption of constant pressure on the conical portion of the body. This 
causes the overestimation of the entropy swallowing distance. For small bluntness, 
the overexpansion can be neglected. Thus, Rotta’s results are valid only for very 
small values of nose bluntness. The distances to attain conical values of surface 
pressure, shock standoff distance, wall heat transfer and skin friction are given in 
Tables 5.3 and 5.4. The effects of cone angle has been investigated for a freestream 
Mach number of 20 and cone angles of 5°, 10° and 20°. Some results for surface 
quantities are presented in Figs. 5.9 - 5.12. An interesting general trend is noticed 
from these figures. With a small increase in cone angle from 5° to 10°, the merging 
distances decrease dramatically. The merging distances decrease even more when 
the cone angle is increased from 10° to 20°. This decrease in merging distance 
can be attributed to the increase in shock angle with increasing cone angle thereby 
reducing the entropy layer thickness. 
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Table 5.2 Entropy swallowing distance (x/iZ n ) for 5° cone at = 10 


Rn 

Reynolds No. 

Present 

Rotta 

0.0025 

2891 

139 

224 


14454 

206 

375 

0.025 

22908 

243 

470 


57816 

286 

590 

0.05 

57816 

288 

590 
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Table 5.3 Distance (x/i?*) up to which bluntness effect persist for M m = 10 


Quantity 

Rn,m 

* e = 5» 

Ogive 
(l = 1-4) 

P* 

0.025 

120 

47 


0.05 

160 

48 



0.025 

209 

63 

S.H 





0.05 

218 

71 



0.025 

233 

127 

Cf 

0.05 

267 

154 



0.025 

120 

44 

Qw 





0.05 

160 

48 
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Table 5.4 Distance (x/-ft„) up to which bluntness effect persist for M - 20 


Quantity 

Rn,m 

6 C = 5° 

O 

o 

II 

o 

0 e = 20° 

Ogive 

(1 = 1.4) 

Ogive 

(EQUIL) 


0.025 

123 

38 

9 

48 

38 

P* 

0.05 

146 

40 

15 

52 

40 



0.025 

172 

43 

10 

57 

44 

$8k 


0.05 

182 

45 

11 

67 

50 



0.025 

165 

80 

25 

87 

65 

Ct 


0.05 

198 

108 

35 

121 

119 



0.025 

126 

37 

10 

46 

36 

Qu, 


0.05 

135 

40 

18 

50 

37 
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The variation in the wall pressure with axial distance for various cone 
angles is shown in Fig. 5.9. The wall pressure increases with the increasing cone 
angle due to increased shock strength. Also the overexpansion/recompression in- 
creases with the cone angle while the physical distance over which this takes place 
decreases. Thus a blunt cone with a larger cone angle will experience a constant 
wall pressure over most of its surface as compared to a cone with a smaller angle. 
Figure 5.10 shows the variation of shock standoff distance with the axial distance 
for 20° cone. The shock standoff distance decreases with the increasing cone angle 
over forward portion. Also the distance it takes to attain conical value decreases 
with the increasing cone angle. Figure 5.11 shows the variation of wall heat transfer 
for various cone angles. It is seen that the heat transfer on the wall increases with 
increasing cone angle while merging distances decrease. Thus larger angle cones 
will experience almost uniform heat transfer over most of the afterbody. Similar 
behavior can be seen for the skin friction in Fig. 5.12. These results indicate that 
as the cone angle increases the extent of effects of nose bluntness decreases. 

The variation of the wall quantities with the Reynolds number is shown 
in Figs. 5.13 - 5.16. In order to change the Reynolds number, the nose bluntness is 
kept constant as 0.025m, while the freestream pressure is changed. The Reynolds 
number is based on the nose radius. The variation of wall pressure and shock 
standoff distance is shown in Figs. 5.13 and 5.14. As expected, these quantities 
were found to be independent of the Reynolds number. The variation of the wall 
heat transfer and the skin friction is shown in Figs. 5.15 and 5.16. A large variation 
can be seen with increasing Reynolds number. With increasing Reynolds number, 
the skin friction and the heat transfer decrease. This is in agreement with the 
laminar theory which states that Cf and C\ are inversely proportional to \/Re. 

Figures 5.17 - 5.23 show the effect of nose bluntness on ogive. Calculations have 
been made for three nose bluntnesses and two freestream conditions. Qualitatively, 
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these results are very similar to those for the cone, although the merging distances 
have decreased. These results are also given in Tables 5.3 and 5.4. The variation 
of wall pressure with axial distance is shown in Fig. 5.17. Here, the sharp ogive 
values are reached earlier than sharp cone values. Unlike cone, the pressure in this 
case does not attain a constant value but decreases gradually with a constant slope. 
The bluntness effects vanish earlier for Mach 20 than for Mach 10 and the region 
influenced by the nose has a lower pressure than that on a sharp body. Figure 5.18 
shows the variation of the shock standoff distance. As expected, the shock is closer 
to the body for Mach 20 than for Mach 10. The skin friction and wall heat transfer 
are shown in Figs. 5.19 and 5.20 respectively. In this case also the skin friction is 
found to be most sensitive to the nose bluntness. It is the last quantity to attain 
the sharp ogive values. 

Figure 5.21 compares the variation of the wall pressure at M = 20 for 5° 
and 10° cone with the ogive for two values of nose bluntnesses. The wall pressure 
for the cone attain a constant value irrespective of the cone angle and bluntness at 
certain distance downstream. But for the case of ogive, the pressure continuously 
decreases, i.e., it does not attain a constant value as in the case of cone for the 
range of axial distances considered in the present study. 

Figures 5.22 and 5.23 show the variation of the streamwise velocity and 
temperature as a function of y distance for Mach 20 perfect gas flow over an ogive 
with various nose bluntnesses. It should be noted that since the shock fitting proce- 
dure is used, the normal profiles are plotted within the shock layer only. The velocity 
and temperature have been nondimensionalized with respect to their freestream val- 
ues. The variation of streamwise velocity is shown in Fig. 5.22 for x = 0.5, 1.0, and 
5.0 m. There is a significant difference between the profiles near the nose; differ- 
ences tend to diminish with increasing axial distance. The figure shows that the 
thickness of the boundary layer increases with increasing bluntness and the shape of 
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the profile is changed near the outer edge of the boundary layer. It is this difference 
which causes the wall values to be so different from the sharp cone values. The 
same type of effects can be seen in Fig. 5.23 where temperature profiles between 
the shock and the body are shown. Here again the differences decrease with the 
axial distance. 

The results for perfect and equilibrium gas flows are shown in Figs. 5.24 
- 5.28. The calculations were performed for Mach 10 and 20, but the results for 
Mach 10 were found to be identical for perfect and equilibrium gas, so only Mach 
20 results are presented. 

The variation of wall pressure, shock standoff distance, heat transfer, and 
skin friction at Mach 20 for ogive of 0.05m nose radius are shown in Figs. 5.24 - 
5.27 for perfect and equilibrium gas flows. The equilibrium gas results do not show 
any dramatic differences from the perfect gas results although, as expected, there 
are some minor differences. The difference in two solutions tends to increase with 
increasing Mach number and bluntness. The most severe case is shown here. These 
figures show expected effects of equilibrium gas chemistry on the flow properties. 

Figure 5.24 shows that the surface pressure variation is almost identical 
for the two cases, except in the recompression zone, where the equilibrium pressure 
is slightly higher than the perfect gas value due to higher density. The variation of 
shock standoff distance is shown in Fig. 5.25; the shock moves closer to the body 
for equilibri um flow calculations. The shock curvature is slightly reduced, thus the 
entropy generation is reduced too. This in turn, decreases the merging distances 
as can be seen in Table 5.4. The decrease in the shock standoff distance can be 
attributed to the increase in density of the gas behind the shock. The variation of 
the skin friction and the heat transfer is shown in Figs. 5.26 and 5.27 respectively. 
The skin friction is reduced slightly in the overexpansion zone but reverses its trend 
in the recompression zone before settling down to the perfect gas value where the 
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Fig. 5.22 Comparison of velocity profiles at x = 0.5, 1.0, and 5.0 m for ogive at 
Mach 20. 
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real gas effects are practically zero. In general, the real gas effects do not affect the 
wall quantities to any large extent for the present conditions. The only noticeable 
difference is for the case of shock standoff distance. 

Figure 5.28 shows the effect of equilibrium gas chemistry on the temper- 
ature profile. It should be noted that even though the maximum temperature in 
the shock layer decreases for the equilibrium gas as compared to the perfect gas, 
the gradients near the wall do not change significantly and therefore the wall heat 
transfer is not affected much. The merging distances for these cases are given in 
Table 5.4. Due to the shock moving closer to the body, a small decrease in merging 
distances is noted for equilibrium flows. 
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Fig. 5.28 Comparison of temperature profiles for equilibrium and perfect gas at x 
= 0.5, 1.0, and 5.U m tor ogive at Mach 20. 



101 


5.2 Forebody and the Cowl Shock Interaction 

In this section the results of the forebody and the cowl shock inter- 
action are discussed. A two-dimensional model is used to study the interaction. 
The computational domain is shown in Fig. 5.29; the solution over the whole do- 
main is obtained in two steps. First, the cowl forebody solution is obtained using 
the Navier-Stokes equations. This solution is then used as the starting solution 
for the parabolized Navier-Stokes equations in the second step to calculate the ex- 
tent of downstream effects of the shock-shock interaction into the inlet. Two sets 
of freestream conditions at Mach 8.03 and 5.94 (Table 5.5) are considered due to 
availability of the experimental data on the cowl forebody. However, no experimen- 
tal data is available for the cowl afterbody. The results for Mach 8.03 are discussed 
first, followed by Mach 5.94 results. For all cases, the undisturbed blunt body flow 
was calculated first and this solution was then used as initial condition for shock 
impingement calculations. First the results are presented for unadapted grid and 
then with the adapted grid to show the importance of a suitable grid in obtaining 
improved solution. 

Figures 5.30 - 5.33 show the temperature, Mach number, pressure, and 
density contours respectively for the cowl forebody at freestream Mach number 
8.03 with unadapted grid shown previously in Fig. 4.10. The surface pressure 
and the heat transfer results are shown in Figs. 5.34 and 5.35 respectively. The 
improved results with adapted grid are shown in Figs. 5.36 - 5.43. The temperature 
and Mach number contours for the cowl afterbody are shown in Figs. 5.43 - 5.46. 
The variation of shock standoff and surface pressure on the cowl afterbody are 
shown in Figs. 5.47 and 5.48. The leeward and windward side computations were 
made independent of each other due to the two-dimensional nature of the flow. 
Similar calculations for Mach 5.94 are shown in Figs. 5.49 - 5.57. For this case the 
calculations were made for adapted grid only. 
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The temperature, Mach number, pressure, and density contours are shown 
in Figs. 5.30 - 5.33, respectively. The distorted bow shock is clearly evident in the 
figures. The bow shock has moved toward the body on the windward side and away 
from it on the leeward side. Here, the windward side is defined as the upper side 
where the flow passes through the impinging shock before encountering the body. 
The interaction of the bow shock and the impinging shock produces a supersonic 
jet surrounded by subsonic flows. The jet terminates with jet shocks and impinges 
on the body producing a local zone of very high pressure and heating rates. The 
stagnation point moves toward the windward side. The location of the stagnation 
point depends upon the strength and orientation of the impinging shock. Due to the 
coarseness of the grid away from the wall, the shocks and shear layers are smeared 
over several grid points. 

The variation of wall pressure along the cowl forebody surface is shown 
in Fig. 5.34. The pressure is nondimensionalized by the stagnation point pressure 
for unimpinged blunt body flow. In order to properly visualize the effect of shock- 
shock interaction, the surface pressure for unimpinged blunt body flow (for the same 
freestream conditions) is also shown in the figure. The pressure on the windward 
side increases considerably with a localized zone of high pressure and falls below the 
surface pressure for unimpinged case on the leeward side. The peak value of pressure 
is about nine times the stagnation point value. The results compare very well with 
the experimental data of Wieting [34] and numerical calculations of Klopfer [44]. 

Figure 5.35 shows the heat transfer along the wall on the cowl forebody. 
The heating rates are nondimensionalized with respect to the stagnation point heat- 
ing for unimpinged case. It shows similar behavior as the surface pressure, i.e., an 
increased heating on the windward side and decreased heating on the leeward side. 
A localized zone of intense heating is observed on the windward side. The results 
compare fairly well with the results of [44] but show poor comparison with the 
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Table 5.5 Freestream and stagnation point conditions for the cowl. 



8.03 387500 2.382 0.0762 12.5 83.5 6.80 E + 5 

5.94 186000 6.857 0.0254 15.0 45.8 2.16 E + 4 
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Fig. 5.30 Temperature contours for the cowl forebody, Moo — 8 03. 
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Fig. 5.32 Pressure number contours for the cowl forebody, Moo — 8.03. 
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Fig. 5.33 Density contours for the cowl forebody, A/oo = 8.03. 
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experimental data. The trend is predicted but the magnitude is relatively lower, 
specially near the jet impingement location. It should be noted that for this case, the 
solution was found to be marginally unsteady as the residual did not go down more 
than three orders of magnitude. Also the jet oscillated slightly, thereby changing 
the location of the peak pressure and heating. Similar behavior was noted for the 
Type IV interaction in [44] and [45]. Due to poor comparison of heating rates and 
smearing of important flow features, it was decided to use adaptive grids thereby 
placing the grid points where they are needed the most. Body curvature, pressure 
and density were chosen as weight functions to adapt the grid to the solution. The 
body curvature clusters the grid points very close to the body while pressure and 
density attract the grid points near shocks and shear layers. Thirty percent of the 
points were allocated for adaptation by body curvature, forty percent by pressure 
and density and the rest of the points were used for creating uniformity of the grid 
so that the grid is not too coarse in any section. The adapted grid is shown in Fig. 
4.10. 

The temperature, Mach number, pressure and density contours are shown 
in Figs. 5.36 - 5.39, respectively for the adapted grid. Since the grid is fine near the 
shocks and shear layer, the flow features are now captured very well. The location 
of the stagnation point, jet and shear layer originating from the shock intersection 
can be seen clearly. The shear layer originating from the stagnation zone is much 
thicker on the leeward side as compared to the windward side. 

Figure 5.40 shows the velocity vectors for the cowl forebody flow. Even 
though the flow is at zero angle of attack, its direction is changed as it passes 
through the impinging shock. The shear layer on the leeward side is seen clearly 

originating from the stagnation point. 

The variations of pressure and heating rate along the cowl forebody surface 
are shown in Figs. 5.41 and 5.42, respectively. The general features of the pressure 
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Fig. 5.38 Pressure contours for the cowl forebody for Moo = 8.03 with adapted grid. 
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distribution are unchanged, except for the peak pressure which is now the same as 
obtained in [44]. It should be noted that for this case the location of the impinging 
shock was made to coincide with the experimental location by matching the peak 
pressure location. The heating rates show a remarkable improvement over the 
previous calculations. The results compare favorably with the experiment. The 
discrepancy in the peak values is probably due to the unsteadiness in the flow and/or 
the turbulent nature of the jet as indicated recently by Moon et al. [46]. The rate 
of convergence for this case was found to be slower than that for the unadapted 
grid case. The cowl of the scramjet engine is modeled by a 5° wedge. The cowl 
afterbody flow is computed by solving the parabolized Navier-Stokes equations using 
the van Leer scheme described previously. The initial plane solution, to initiate the 
marching procedure, is provided by the cowl forebody solution. The objective of 
carrying out the cowl afterbody calculation is to determine how the flow entering 
the engine inlet is affected by the shock-shock interactions. 

Figures 5.43 and 5.44 show the temperature contours on the cowl afterbody 
for the windward and the leeward sides, respectively. The flow field on the windward 
side is quite different as compared to the leeward side. The flow on the windward 
side shows typical two-dimensional wedge flow features, such as rays of the constant 
Mach number and the shock attaining conical angle slightly downstream of the nose. 
The leeward side shows thick shear layer interacting with the boundary layer; this 
interaction again produces a weak shear layer whose strength diminishes rapidly 
with the axial distance. Similar behavior can also be seen in the Mach contours for 

the cowl afterbody in Figs. 5.45 and 5.46. 

Figure 5.47 shows the variation of shock standoff distance along the cowl 
afterbody on leeward and windward sides for the impinged case. These are com- 
pared with the results for the unimpinged case. The shock moves closer to the 
body on the windward side and away from it on the leeward side as compared to 
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Fig. 5.42 Variation of surface heat transfer for the cowl forebody, Moo = 8.03 with 
adapted grid. 


120 


the unimpinged case. The calculations were carried out up to 40 nose radii down- 
stream. The angle at which the shock enters the inlet is significantly changed due 
to the shock-shock interaction which can affect the inlet performance. 

The variation of pressure on the cowl afterbody is shown in Fig. 5.48. 
The pressure on the windward side increases on the forward portion of the cowl 
afterbody; it reaches a peak value and then settles down to a constant value as the 
flow recompresses back after the expansion on the cowl forebody. On the leeward 
side also, the pressure increases and then gradually attains approximately the same 
value as for the unimpinged case. The rate of pressure increase on the forward 
portion of the leeward side is found to be slower than on the windward side. 

Similar calculations were carried out for Mach 5.94 freestream conditions 
and the results are shown in Figs. 5.49 - 5.58. Qualitatively these results are very 
similar to the Mach 8.03 conditions. The grid for this case is generated in the same 
way as for the previous case. The temperature, Mach number, pressure and density 
contours are shown in Figs. 5.49 - 5.52 for the cowl forebody. As noted previously 
the shock moves toward the body on windward side and away from it on leeward 
side. The effects of jet are clearly visible in the figure. The variations of surface 
pressure and heat transfer are shown in Figs. 5.53 and 5.54, respectively. The 
numerical results are compared with the experimental data reported by Tannehill 
et al. [42] from an unpublished experiment by J. W. Keyes of NASA Langley. The 
computed values of surface pressure agree very well with the experimental data 
although the peak pressure is slightly lower than the experimental peak. Although 
not shown here, the peak value of the pressure matches with the numerical calcu- 
lations of Tannehill et al. [42]. The variation of heat transfer on the cowl forebody 
is shown in Fig. 5.54. For heating rates, the experimental data was reported only 
at one point. The peak value of the heating rate is 6.4 times the stagnation point 
value, which is within the uncertainty range of the experiment. 
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The results for the cowl afterbody are shown in Figs. 5.55 - 5.58. The 
Mach number contours for the windward and the leeward sides on the afterbody are 
shown in 5.55 and 5.56. The shear layer on the leeward side is thinner as compared 
to the Mach 8.03 due to weaker strength of the interacting shocks. The variation 
of shock standoff distance for the cowl afterbody is shown in Fig. 5.57 and the 
variation of the wall pressure for the cowl afterbody is shown in Fig. 5.58. The 
results show the same trends as for the previous case. 









Fig. 5.47 Variation of shock standoff distance with axial distance for the cowl 
afterbody, A/qo — 8.03 



Fig. 5.48 Variation of surface pressure with axial distance for the cowl afterbody 
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Fig. 5.53 Variation of surface pressure for the cowl forebody, Moo — 5.94 
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Fig. 5.54 Variation of surface heat transfer for the cowl forebody, = 5.94. 
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Fig. 5.57 Variation of shock standoff distance with axial distance for the cowl 
afterbody, Moo = 5.94 
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Fig. 5.58 Variation of surface pressure with social distance for the cowl afterbody, 
Moo = 5.94 
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5.3 Forebody Shock and the Inlet Sidewall Shock 

Interaction 

In this section the interaction of the forebody shock with the inlet sidewall 
shock is discussed. As noted previously, the flow is fully three-dimensional in nature 
for this case due to the orientation of the impinging forebody shock as well as the 
sweep of the inlet. Due to sweep of the inlet sidewall, the Type V interaction occurs 
and is considered in this section. Results have been obtained for a 25° swept sidewall 
and 10° impinging shock generator angle. The freestream conditions used for the 
present calculations are given in Table 5.6. For three-dimensional calculations also, 
the flow field is initialized by the undisturbed flow field solution. First the results 
are presented for the undisturbed flow and then with the impinging shock to show 
the effect of impingement. 

Figure 5.59 shows the pressure contours for the unimpinged case. The 
stagnation plane ( 0 = 0) contours are shown in Fig. 5.59(a). The bow shock can 
be seen clearly and is swept back as it follows the body contours. The flow in each 
cross-sectional plane is identical, so only one such plane is shown in Fig. 5.59(b). 
The figure shows a typical blunt body solution. At the stagnation point, the pressure 
has its peak value and the shock standoff distance is minimum. Figure 5.60 shows 
the temperature, Mach number and the pressure contours in the stagnation plane 
with shock impingement. The temperature contours are shown in Fig. 5.60(a). 
Due to stretching, the grid away from the surface is coarse; consequently, the shock 
and the other flow features are smeared over several grid points. The impinging 
shock is clearly visible in the figure. The bow shock moves away from its unimpinged 
position and a shear layer is produced . The shear layer interacts with the boundary 
layer. Also, a transmitted shock is produced which strikes the body causing a jump 
in pressure and heating rates. These are typical features of the Type V interaction. 
There are temperature and density gradients across the shear layer. The Mach 
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Table 5.6 Freestream conditions for inlet. 


Moo 

Re D 

T 

± IU 

Too 

Rn 

6 

5.94 

180000 

394 k 

59.6 k 

0.0125 m 

10.0 
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number contours are shown in Fig. 5.60(b); the flow is fully subsonic behind the 
bow shock and the shear layer can be seen. Figure 5.60(c) shows the pressure 
contours in the stagnation plane. Since the pressure is constant across the shear 
layer, it is not visible in this figure. Some expansion waves can be seen emanating 
from the intersection point. These expansion waves reduce the surface pressure as 
the flow expands. 

More insight into the numerical solution can be gained by examining the 
plots for z = constant planes. Figure 5.61 shows a series of contour plots for in- 
creasing values of z. The movement of the impingement point (as it moves away 
from stagnation point) can be seen clearly; also, the shear layer as it develops is 
very well captured. 

Figure 5.62 shows a cylindrical fin coated with temperature sensitive paint 
and the corresponding schlieren photograph locating regions of high pressure and 
heat transfer for Type V interference at Mach 4.6 and is taken from [31]. It should 
be noted that the freestream conditions for this case are slightly different but it 
shows a typical Type V interaction. The numerical solutions are shown in Figs. 
5.63 and 5.64. The surface pressure contours are shown in Fig. 5.63 and the 
stagnation plane pressure contours in Fig. 5.64. The numerical solution captures 
the first two (A and B) peaks in pressure and heat transfer. The third peak (C) is 
caused due to separation near the nose which is absent in the numerical solution due 
to assumption of infinite length of the body. The behavior of bow shock near the 
impingement point is very well captured by the numerical scheme. The comparison 
of Figs. 5.63(b) and 5.64 shows how well the physical phenomenon is capture by 
the numerical scheme. 

Figures 5.65 - 5.70 show the variation of pressure and heat transfer along 
the surface. They are nondimensionalized with respect to the unimpinged stag- 
nation line values. The unimpinged stagnation line pressure is obtained by the 
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Fig. 5.62 Local regions of high heating rates and corresponding schlieren photo- 
graph. 
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Fig. 5.63 Surface pressure contours for Type V interaction. 
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Fig. 5.64 Pressure contours in the stagnation plane(0 — 0). 
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Rayleigh pitot formula as 


Po _ , m ?)h 6 

5 j 7Mf — 1 


■)* 


(5.1) 


Where Mi is the normal Mach number. The stagnation line heat transfer 
is obtained by the following expression developed by Beckwith and Gallagher [77] 

, T aw — T w 2i?ep )00 Mox i fZ TooPo . p 0 >-|i .2) 

q 0 = k <x -—=—( L n K ; 


where 

r aiu = r e (i + r — Af*) 

here r is the recovery factor and is taken as 0.85 and subscript e denotes the bound- 
ary layer edge conditions. 

Figure 5.65 shows the comparison of stagnation line pressure with the 
experimental data of Keyes and Hains [32] and numerical calculations of Holst et 
al. [47]. In the experimental set up, the interaction point was only 3 cm downstream 
from the end of the cylinder (z = 0) and, therefore, there was some relieving effect 
(Fig. 5.62). But in the present calculations, the body is assumed as infinite in 
length and hence no relieving effect is allowed. Due to this discrepancy in the 
model used, the results do not match near the end. However, the results compare 
fairly well away from the end point. The peak value of the pressure is caused by 
the interaction of the transmitted shock with the boundary layer and it is very well 
captured. The flow overexpands and then recompresses back to the unimpinged 
value. The peak value of the pressure is about 2.2 times the unimpinged stagnation 
line pressure. The variation of the stagnation line heat transfer is shown in Fig. 
5.66. In this case also, the comparison is poor near the end for the same reasons as 
explained earlier. But away from the end, results compare fairly well. As expected, 
the heat transfer follows the same general trend as the pressure. The peak value 
of heat transfer is about three times the unimpinged stagnation line heat transfer. 
The experimental data is available only for the stagnation plane. 
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The results from the present calculations show a remarkable improvement 
over Holst’s calculations due to various factors. The flow field is much better re- 
solved due to a finer grid, no artificial thickening of the boundary layer has been 
attempted and the boundary condition at the end( z=0 ) are physically correct. The 
present calculations employ zero gradient conditions at the end which are physically 
realistic. In the previous calculations, the flow was maintained at constant condi- 
tions corresponding to the unimpinged flow solution, which is incorrect for the 
impinging shock case. 

The variation of pressure and heat transfer at the surface is shown in 
Figs. 5.67 and 5.68 for various values of y. The stagnation line (y=0) has the 
maximum amplification of pressure and heating rates as the bow shock has its 
maximum strength along this line. Away from the stagnation line, strength of bow 
shock decreases. Consequently, the interaction becomes weaker and weaker, thereby 
causing less amplification of the wall quantities. 

The circumferential variations of the wall pressure and the heat transfer 
are shown in Figs. 5.69 and 5.70 for various values of z. The impinging shock 
intersects at the normal portion of the bow shock at z = 2.4 cm, and therefore, the 
peak values of pressure and heating occur in this circumferential plane. At higher 
values of z, the impinging shock intersects away from the normal portion of the bow 
shock, thus reducing the amplification. 
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Fig. 5.67 Variation of the surface pressure in various y - constant planes 



Fig. 5.68 Variation of the surface heat tra 




Chapter 6 

CONCLUSIONS 


This study presents a numerical investigation of the effects of blunt lead- 
ing edges on the viscous flow field around a hypersonic vehicle such as the proposed 
National Aero-Space Plane. Attention is focused on two specific regions of the flow 
field. In the first region, effects of nose bluntness on the forebody flow field are 
investigated. The second region of the flow considered is around the leading edges 
of the scramjet inlet where the forebody shock interacts with the shock produced 
by the blunt leading edges of the inlet compression surfaces. 

An accurate solution procedure is developed to determine the downstream 
influence of forebody bluntness and the conditions, under which it can be modeled 
by a relatively simple pointed nose body of revolution, are determined. A wide 
range of physical and geometrical conditions are considered. The geometrical mod- 
els used are spherically blunted cones of 5°, 10° , and 20° angles and ogives with 
various degree of bluntness. Several approximate methods as well as combination 
of Navier-Stokes and parabolized Navier-Stokes equations are used to accurately 
compute the flow field. The high temperature effects are incorporated by using 
equilibrium gas model. The influence of entropy layer generated by the shock curva- 
ture on the flow field are also addressed. The results show that the effects of leading 
edge bluntness persist several hundred nose radii downstream depending upon the 
bluntness, freestream conditions and geometry. It is found that the entropy layer 
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generated by the shock curvature affects the development of the boundary layer far 
downstream of the nose. These effects decrease with increasing Mach number. Also 
the bluntness effects decrease with increasing cone angle. The results show that the 
wall quantities are not affected much by the inclusion of high temperature effects 
through equilibrium chemistry for the conditions considered in the study. 

The interaction of the forebody shock with the blunt cowl shock and with 
the shock produced by the blunt leading edge shock of the inlet is also investigated. 
The calculation of flow field and aerothermal loads under such conditions is quite 
difficult due to complexity of the flow field and the presence of high gradients. A 
modern computational algorithm developed by van Leer which has the advantage 
of better modeling of physics of the problem has been used. The impinging fore- 
body shock is treated as a sharp discontinuity across which the exact shock jump 
conditions (Rankine-Hugoniot) are applied. The compressible Navier-Stokes equa- 
tions are solved using a shock capturing algorithm. It does not require any prior 
information about the shock impingement flow field to be computed as is the case 
with some of the previous semi-empirical approaches. The use of solution adaptive 
grid technique to properly resolve the flow gradient is made and the importance 
of suitable grid in the numerical calculations is shown. The numerical computa- 
tions are primarily done for Type IV supersonic jet interaction since it represents 
the most severe case in terms of aerothermal loads. A preliminary investigation 
of the Type V interaction which occurs on the leading edge of the compression 
side walls of the inlet is also made. The flow field for this interaction is found to 
be fully three-dimensional unlike the cowl plate interaction. Results of the study 
show the presence of local regions of very high surface pressure and intense heating. 
For Type IV interaction, the peak amplification of pressure is found to be nine 
times and heating amplification to be eight times the stagnation point value of the 
unimpinged case. The degree of amplification is found to increase with the increas- 
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ing Mach number. For Type V interaction, the amplification factors are found to 
be much lower, the pressure is amplified by a factor of 2.25 and heat transfer three 
times the stagnation values for unimpinged case. The results compare well with 
the experimental data. The study demonstrates the capability of CFD to simulate 

extremely complex viscous hypersonic flows. 

Based on the present study, several recommendations concerning the ex- 
tensions of this work are suggested. 

(1) The realistic forebody geometry should be used for calculating bluntness effect. 

(2) The flow in the present investigation is assumed to be laminar, but in actuality 
the flow would be turbulent, so an appropriate turbulence model should be incor- 
porated into the governing equations. 

(3) Due to high Mach number flows considered, it would be desirable to examine 
the non-equilibrium high temperature effects. 

(4) The oscillatory behavior of the supersonic jet in Type IV interaction needs to 
be further investigated. 

(5) Since the inlet will have more than one shock impinging on it, the effects of 
multiple shock wave interaction needs to be determined. 

(6) The effects of oscillations in the impinging shock should be investigated. 

(7) The use of adaptive grids for three-dimensional calculations to properly resolve 

the flow should be made. 

(8) For Type V interaction considered, the end effects should be accounted for to 
have a better comparison with the experiment by using zonal techniques. 

(9) The effects of shock wave interference in the presence of active cooling should 
be determined because to relive the thermal stresses, some form of active cooling 

would be used. 

(10) Various shock interference pattern for three-dimensional calculations should be 
computed by varying the sweep angle. 
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APPENDICES 


APPENDIX A 


EQUILIBRIUM CHEMISTRY 


A fully vectorized code EQUIL as described in [56] has been used in 
the present analysis to include the real gas effects. It is based on free energy 
minimization technique [78] in which the method of steepest descent is used. In this 
method, with the known values of local temperature, pressure and elemental mass 
fractions of the gas, the mole numbers of species are calculated which minimizes 
the free energy. The free energy of a mixture of n species containing x, moles of the 
ith species is expressed as 

F{X) = £ /• ( A>1 ) 

»=1 

where 


X — (Xi,X2j^3 x n) 


(A. 2) 


ft — x i 


Ci + In— 
x J 




=x> 


(A. 3) 
(A. 4) 
(A. 5) 


Equation (A.l) is minimized, subject to mass balance constrain. An initial 
assumption is made on the mole numbers of various species, and then an iterative 
procedure is followed to find the set of mole numbers which satisfy both the con- 
straints. Once the mole numbers are known, the molecular weight and mass fraction 
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of each species are calculated. The specific heat, free energy, viscosity and thermal 
conductivity are calculated using polynomial expressions given below. 


^ = a x + a 2 T + o 3 T 2 + a 4 T 3 + a 5 r 4 

H T T 2 T 3 T 4 1 

— = ^ + 02 ^ + a 3 — + a 4 — + a 5 — + a*- 

F T T 2 T s T * 

— = ai (l - InT ) - a 2 - - a 3 — ~ a 4 — - a 5 — + *6 



H = 6i + b 2 T + b 3 T 2 


fc = c x + c 2 T 


(A. 6) 
(A. 7) 

(A. 8) 
(A. 9) 
(A. 10) 


The values of constants ai,a 2 , ,ar,bi, ,b 3 ,c l and c 2 are given in [56] 


APPENDIX B 
ADAPTIVE GRID 


An adaptive grid system was employed to improve the quality of the so- 
lution without increasing the number of grid points. The accuracy of the numerical 
solution depends on the fineness of the grid. The presence of large gradients causes 
the error to be large in numerical approximation of the derivatives. Therefore, it 
is desirable to have fine grid spacing near the shocks, shear layers, and boundary 
layer. The method developed by Abolhassani et al. [76] has been used. It is a very 
general method with capability to adapt the grids with various variables such as 
pressure, Mach number, shear stress etc and is based on variational approach. It 
is an algebraic method and is formulated in such a way that there is no need for 
matrix inversion, which makes it computationally very efficient. The grids were 
adapted by using Euler Lagrange equations 

X^w = constant (®* 


This equation can be written in normalized form as 


t{») 


Jo 

fo w(t)dt 


(B.2) 


where 3 is the arc length. The weight function is used to reduce the grid point 
spacing where w is large. For grid adaptation with several variables, the weight 
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function is expressed as 


u> = i+£>/. < B - 3 ) 

1=1 

where n is the number of variables, 6, are the constants and /< are the variables ( or 
their derivatives ). The constant 1 is for uniformity. The values of 6, are calculated 
based on percentage of grid points allocated to each variable. The percentage of 
grid points for each /, can be written as 


bj Fj ( S m ax ) 

Rj= S max + Z: = 1 b t F t (S maz ) 


-;j = 1,2,3, 4. 


(B.4) 


where 

Jim = [ m<u 

The integral is approximated by a trapezoidal rule. After some algebraic manipu- 
lations, Eqs. (B.2) and (B.4) can be written as 


f(s) = 




+ Eft 


JIM 


»=l 


^ Fi{S ( max) 


(B.5) 


In order to avoid extremely small grid spacing near the shocks, the weight functions 
f i are multiplied by the factor 


1 - (B. 6) 

where A S„i„ is the allowable minimum spacing. The function varies from 0 to 1 
and is proportional to the spacing. 



APPENDIX C 


RADIATIVE INTERACTIONS 
IN LAMINAR FLOWS 


A preliminary study was conducted to include the radiative formulation 
in the general governing equations and provide the step by step analysis and solu- 
tion procedure for realistic problems. The objective of this study was to prepare 
the ground work for future studies of inlet flows involving radiative interactions. 
The original plan was to compute the inlet flow in combination with shock-shock 
interaction as well as radiative heat transfer. But due to time limitation this part of 
study could not be completed. However, some very important results were obtained 
during the course of this study. A brief introduction to the problem is given here, 

the details are available in [79] and [80]. 

The problem has been formulated in such a way that any nongray gas 

model for participating species can be used. The governing equations involving ra- 
diative interactions, in general, are integro-differential equations. In order to fully 
understand the radiative part, it was decided to use a very simple fluid dynamic 
model along with sophisticated nongray gas models. The laminar flow of participat- 
ing species between two parallel plates was considered. The Tien and Lowder [81] 
wide-band correlation model was used to account for the spectral properties. The 

walls are assumed as black. 
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Specific results were obtained for energy transfer by pure radiation, ra- 
diation and conduction, and finally by radiation, conduction and convection. The 
participating species considered were CO, OH, C02, H20 and their mixtures. It 
was noted that the extent of radiative interaction is dependent on the nature of 
participating species and parameter such as Tl, T2/T1, P and PL. Here T1 and 
T2 are lower and upper wall temperatures, P is the pressure and L is the plate 
spacing. Radiative ability of a gas was found to increase with increasing pressure, 
temperature and plate spacing. For most conditions, H20 was found to be a highly 

participating species (as compared to CO, OH and C02) 

The analysis, computational procedure, and results presented in [79] and [80] 
provide essential information on how to account for radiative interaction in internal 
flows. The models developed in the references can be used to investigate radiative 
interaction in complex flows as long as the system temperature remains within the 


range 300 - 5,000 K. 



